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1 Introduction 

In this paper we address the question of stability of Minkowski space-time for the system of the 
Einstein-scalar field equation^] 

(1-1) — ^9fiuR = T^ v . 

The equations connect the gravitational tensor G^ u = R^ u — \g^v given in terms of the Ricci R^y and 
scalar R = g^Rnv curvatures of an unknown Lorentzian metric g^y and the energy- momentum tensor 
T^y of a matter field ip: 



X2) ^=8, 



^ip d v ip - ^g^y (g a ^d a ip dptp) 



The Bianchi identities 



D^G^y = 



imply that the scalar field tp satisfies the the covariant wave equation 

□ fl V = Tfj M Jg^^\fetg~\dyA = 

The set (m,R 3+1 ,0): standard Minkowski metric g = m = —dt 2 + X]i=i(^ x *) 2 on I^ 3+1 an d vanishing 
scalar field ifi = describes the Minkowski space-time solution of the system (jl.ip . 

The problem of stability of Minkowski space appears in the Cauchy formulation of the Einstein 
equations in which given a 3-d manifold T,q with a Riemannian metric go, a symmetric 2-tensor ko 
and the initial data (tpo, tp\) for the scalar field, one needs to find a 4-d manifold Ai, with a Lorantzian 
metric g and a scalar field ip satisfying the Einstein equations (II. lj) . and an imbedding So C M such 
that go is the restriction of g to S, ko is the second fundamental form o/E and the restriction of tp to 
So gives rise to the data (ipo,ip\). 



*Part of this work was done while H.L. was a Member of the Institute for Advanced Study, Princeton, supported by 
the NSF grant DMS-0111298 to the Institute. H.L. was also partially supported by the NSF Grant DMS-0200226. 

tPart of this work was done while I.R. was a Clay Mathematics Institute Long-Term Prize Fellow. His work was also 
partially supported by the NSF grant DMS-01007791. 

1 We use Greek indices a, /3, fj,, v... — 0, ...,3, the summation convention over repeated indices and the notation d a = 
d/dx a . The symbol D denotes a covariant Levy-Civita derivatives with respect to the metric g 
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The initial value problem is over determined and the data must satisfy the constraint equations: 

Ro - h) koj + kol koj = |VVo| 2 + iV'il 2 , V^oii " v * **>j = V^o fh- 

Here Rq is the scalar curvature of go and V is covariant differentiation with respect to go- 

The seminal result of Choquet-Bruhat [CBlJ followed by the work [CB-GJ showed existence and 
uniqueness (up to a diffeomorphism) of a maximal globally hyperbolic!! smooth space-time arising from 
any set of smooth initial data. The work of Choquet-Bruhat used the diffeomorphism invariance of 
the Einstein equations which allowed her to choose a special harmonic (also referred to as a wave 
coordinate or de Donder) gauge, in which the Einstein equations become a system of quasilinear wave 
equations on the components of the unknown metric g^ v 

(1.3) U g g^ = F tlu {g){dg,dg) + 2d^d^ ) 5^ = 0, where U g = g^d a d 

with F(u)(v,v) depending quadratically on v. Wave coordinates {x' 1 } /J= o,...,3 are required to be so- 
lutions of the wave equations U g x^ = 0, where the geometric wave operator is \3 g = D a D a = 
g a Pd a dp + g a ^T £nd v . The metric g^ u relative to wave coordinates {x^} satisfies the wave coordinate 
condition 

(1-4) 9 aP 9^T^ = g aP dpg afl - \g^d^g a p = 0. 

Under this condition the geometric wave operator \D g is equal to the reduced wave operator U g . The 
use of harmonic gauge goes back to the work of Einstein on post-Newtonian and post-Minkowskian 
expansions. 

In the PDE terminology the result of Choquet-Bruhat corresponds to the local well-posedness of the 
Cauchy problem for the Einstein-vacuum (scalar field) equations with smooth initial data. The term 
local is appropriate in the sense that the result does not guarantee that the constructed space-time is 
causally geodesically complete^ and thus could "terminate" in a singularity. Our experience suggests 
that global results require existence of conserved or more generally monotonic positive quantities. The 
only known such quantity in the asymptotically flat case is the ADM mass, whose positivity was 
established by Schoen-Yau |S-Y] and Witten |Wi| . is highly supercritical relative to the equations 
and thus not sufficient to upgrade a local space-time to a global solution. This leaves the questions 
related to the structure of maximal globally hyperbolic space-times even for generkjfl data firmly in 
the realm of the outstanding Cosmic Censorship Conjectures of Penrose. Given this state of affairs 
the problem of stability of special solutions, most importantly the Minkowski space-time, becomes of 
crucial importance. 

Stability of Minkowski space-time for the Einstein-vacuum (scalar field) equations: 

Show existence of a causally geodesically space-time asymptotically "converging" to the Minkowski 
space-time for an arbitrary set of smooth asymptotically flat initial data (Eoj ffOyj fcoij) with So ~ M 3 , 

(1.5) goij = (1 + — )Sij + o(r- l ~ a ), k 0ij = o(r- 2 ~ a ), r = \x\ -> oo, a > 

2 A space-time is called globally hyperbolic if every inextendable causal curve intersects the initial surface E once and 
only once. A causal curve x(s) is a curve such that g a px a x^ < 0, where x = dx/ds. It is future directed if x° > 0. 
3 i.e. any causal geodesies x(s), g a px a x^ = const < can be extended to infinite parameter value < s < oo. 
4 Note that all the known explicit solutions, with exception of the Minkowski space-time, are in fact incomplete. 
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where {qq — 5) and k$ satisfy global smallness assumptions . The stability of Minkowski space-time for 
the Einstein-scalar field equations: in addition requires a global smallness assumption on the scalar 
field data (tpo,t/ji), which obey the asymptotic expansion 

(1.6) ^ = ( r -l-<*), ^ 1 = ( r -2-a)_ 

A positive parameter M in the asymptotic expansion for the metric go is the ADM mass. 

The stability of Minkowski space for the Einstein-vacuum equations was shown in a remarkable 
work of Christodouolou-Klainerman for strongly asymptotic initial data (the parameter a > 1/2 in the 
asymptotic expansion (|1.5p .) The approach taken in that work viewed the Einstein- vacuum equations 
as a system of equations 

D a W a(3yS = 0, D a * W a0l5 = 

for the Weyl tensor W a /s^s of the metric g a p and used generalized energy inequalities associated with the 
Bel-Robinson energy-momentum tensor, constructed from components of W, and special geometrically 
constructed vector fields, designed to mimic the rotation and the conformal Morawetz vector fields of 
the Minkowski space-time, i.e., "almost conformally Killing" vector fields of the unknown metric g. The 
proof was manifestly invariant, in particular it did not use the wave coordinate gauge. This approach 
was later extended to the Einstein-Maxwell equations by N. Zipser, [Z] . 

Nevertheless the PDE appeal (for the other motivations see the discussion below) of the harmonic 
gauge for the proof of stability of Minkowski space-time lies in the fact that the latter can be simply 
viewed!! as a small data global existence result for the quasilinear system (jl.3p . However, usefulness of 
the harmonic gauge in this context was questioned earlier and it was suspected that wave coordinates 
are "unstable in the large", [CB2| . The conclusion is suggested from the analysis of the iteration 
scheme for the system (|1.3|) : 

9nv = m^u + eg$ + e 2 g$ + 

where g^},gj^} satisfy respectively homogeneous and inhomogeneous wave equations on Minkowski 
background 

(1.7) □<$=(), Ug$=F{jn^){dg$,dg$), where □ = m°^a a dp 

As a solution of the homogeneous 3 + 1-d wave equation with smooth decaying initial data the functions 
gtiv ~ Cet^ 1 in the so called wave zone t ~ r as t — > oo. Integrating the second equation implies 
that gfy ~ Cet^lnt. This suggests that already the asymptotic behavior of the second iterate 
deviates from the one of the free waves in Minkowski space-time and plants the seeds of doubt about 
global stability of such a scheme. We should note that this formal iteration procedure is termed the 
post-Minkowskian expansion and plays an important role in the study of gravitational radiation from 
isolated sources, see e.g. (Blj, [B-DJ. 

To understand some of the difficulties in establishing a small data global existence result for the 
system (|1.3p let us consider a generic quasilinear system of the form 

(1.8) Ufa = b{ kafS d a (j)j dp<f) k + ^2 4^ <Aj d a d/3(t) k + cubic terms 



5 This statement requires additional care since a priori there is no guarantee that obtained "global in time" solution 
<7 M „ defines a causally geodesically complete metric. However, the latter can be established provided one has good control 
on the difference between g M „ and the Minkowski metric L-R2 . 
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The influence of cubic terms is negligible while the quadratic terms are of two types, the semilinear 
terms and the quasilinear terms, each of which present their own problems. The semilinear terms can 
cause blow-up in finite time for smooth arbitrarily small initial data, as was shown by John[JlJ, for 
the equation □ <p = (dt<ft) 2 . D. Christodoulou [Cl| and S. Klainerman |K2j showed global existence for 
systems of the form f|l .8j) if the semilinear terms satisfy the null condition and the quasilinear terms 
are absent. The null condition, first introduced by S. Klainerman in [KlJ, was designed to detect 
systems for which solutions are asymptotically free and decay like solutions of a linear equation. It 
requires special algebraic cancellations in the coefficients b{ kal3 , e.g. □ 4> = (dt4>) 2 — (V^l 2 . However, 
the semilinear terms for the Einstein equations do not satisfy the null condition, see [CB3j. The 
quasilinear terms is another source of trouble. The only non-trivial example of a quasilinear equation 
of the type (|1.8p . for which the small data global existence result holds, is the model equation □</> = 
4>A(j), as shown in |L2j (radial case) and [A3] (general case). However the solutions only decay like 
cf> ~ et~ 1+Ce . This in particular implies that the characteristic surfaces of the associated Lorentzian 
metric g = —dt 2 + (1 + eft) X^=i(^ x *) 2 diverge distance ~ t Ce from the Minkowski cones. 

In our previous work |L-R1| we identified a criteria under which it is more likely that a quasilinear 
system of the form ()2.9|) has global solutions^. We said that a system of the form (|2.9p satisfy the 
weak null condition if the corresponding asymptotic system (c.f. |Hll IH2] ) has global solutiongj We 
showed that the Einstein equations in wave coordinates satisfy the weak null condition. In addition 
there is some additional cancelation for the Einstein equations in wave coordinates that makes it better 
than a general system satisfying the weak null condition. The system decouples to leading order, when 
decomposed relative to the Minkowski null frame. An approximate model that describes the semilinear 
terms has the form 

□02 = (cW>l) 2 , □&=(). 

While every solution of this system is global in time, the system fails to satisfy the classical null 
condition and solutions are not asymptotically free: <p2 ~ hr \t\- The semilinear terms in Einstein's 
equations can be shown to either satisfy the classical null condition or decouple in the above fashion 
when expressed in a null frame. The quasilinear terms also decouple but in a more subtle way. The 
influence of quasilinear terms can be detected via asymptotic behavior of the charachteristic surfaces 
of metric g. It turns out that the main features of the characteristic surfaces at infinity are determined 
by a particular null component of the metric. The asymptotic flatness of the initial data and the wave 
coordinate condition (|1.4|) give good control of this particular component, i.e., ~ M/r, which in turn 
implies that the light cones associated with the metric g diverge only logarithmically ~ Mint from 
the Minkowski cones. 

In |L-R2] we were able to establish the "small data global existence" result for the Einstein- vacuum 
equations in harmonic gauge for special restricted type of initial data. In addition to the standard 
constraint compatibility and smallness conditions the initial data was assumed to coincide with the 
Schwarzschild data 

M 

90ij = (1 + — k 0ij = 

in the complement of the ball of radius one centered at the origin. The existence of such data was 
recently demonstrated in |Co| , |C-D| . The stability of Minkowski space-time for the Einstein- vacuum 

6 At this point, it is unclear whether this criteria is sufficient for establishing a "small data global existence" result for 
a general system of quasilinear hyperbolic equations. 

7 The condition trivially holds for the class of equations satisfying the standard null condition. For more discussion of 
the weak null condition and a general intuition behind the proof see Section 2. 
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equations for such data of course also follows from the work of Christodoulou-Klainerman and yet 
another approach of Friedrich [FY] • This choice of data allowed us to completely ignore the problem of 
the long range effect of the mass and the exterior existence. 

In this paper we prove stability of Minkowski space-time for the Einstein- vacuum and the Einstein- 
scalar field equations in harmonic gauge for general asymptotically flat initial data (any a > in 
p.5jl - (jl.6p ) close to the data for the Minkowski solution. 

The asymptotic behavior of null components of the Riemann curvature tensor R a p^s of metric g- 
the so called "peeling estimates"- was discussed in the works of Bondi, Sachs and Penrose and becomes 
important in the framework of asymptotically simple space-times (roughly speaking, space-times which 
can be conformally compactified) , see also the paper of Christodoulou |C2j for further discussion of 
such space-times. The work of [C-K] provided very precise, although not entirely consistent with 
peeling estimates, analysis of the asymptotic behavior of constructed global solutions. However, global 
solutions obtained by Klainerman-Nicolo |K-Nlj in the problem of exterioi^l stability of Minkowski 
space were shown to possess peeling estimates for special initial data, [K-N2J. 

Our work is less precise about the asymptotic behavior and is focused more on developing a rel- 
atively technically simple approach allowing us to prove stability of Minkowski space in a physically 
interesting wave coordinate gauge, for general asymptotically flat data, and simultaneously treating 
the case of the Einstein equations coupled to a scalar field. 

Theorem 1.1. Let (£, go, ko, V'O) V'l) be initial data for the Einstein- scalar field equations. Assume 
that the initial time slice £ is diffeomorphic to R 3 and admits a global coordinate chart relative to 
which the data is close to the initial data for the Minkowski space-time. More precisely, we assume 
that the data (go,ko,ipo,^i) is smooth asymptotically flat in the sense of (|1.5|l - (|1.6p with mass M and 
impose the following smallness assumption: Let 

go = 5 + hi + hi, where hf jQ = x(r) — <% 
where %(s) G C°° is 1 when s > 3/4 and when s < 1/2. Set 

(1.9) e n (o)= £ iki + ^i/^+t+i/Ivv^ii^ + iki + ov^t+uiv^oIi^ 

0<|/|<7V 

+||(1 + r) 1 /2+7+|/| V V%|| La + ||(1 + r) 1 /2+7+]/| v ^ 1 || L2 
There is a constant eo > such that for all e < Eo and initial data verifying the condition 

(1.10) E N (0) + M<e 

for some 7 > 70(^0) with 7o(eo) ~~ > as eo — > and N > 6, the Einstein- scalar field equations possess a 
future causally geodesically complete solution {g, ip) asymptotically converging to Minkowski space-time. 
More precisely, there exists a global system of coordinates (t, 3 ) in which the energy 

(1.11) 

£ N (t)= Ww^dZ'h^tr^ + Ww^dZ 1 ^,-)^, w l ' 2 = 

\I\<N 

8 Outside of the domain of dependence of a compact set 



(l + |r-f|) 1 /2+7 ) 



r > t 
r < t 
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with Minkowski vector fields Z 6 {d a , x a dp — xpd a ,x a d a } of the solution 

g(t) = m + h°(t) + h\t), h% = x (r/t) X (r)y6 a p. 

obeys the estimate 

(1.12) £ N (t) <C N e(l + tf N£ . 

Moreover, 

(Lis) | 8Z vi + | 8 z^<( c w i+(+rr ;:^ (i+|( - r|) ";;:' r>( - ik*-* 

1,1 11 ' - (CVU + t + rr'+Cwfl + lf-rl)- 1 / 2 , r<(, M ~ 
and 

(i.i4) \ Z ^\ + izvi < |^ £(1 + 1 + rr ;i"r (1 + " - r|) " 7 ' r **> 111 - N ~ 2 

K ' ' 11 ' - \c^e(l + t + r)- 1+2 ^ £ , r<t, |/|<7V-3 

For |/| = we prove the stronger bound Cet" 1 In \t\ and in fact show that all but one component 
of h 1 , expressed relative to a null frame, can be bounded by Cet^ 1 . This more precise information is 
needed in the proof. 

As we will show in a future paper, the decay estimates above are sufficient to prove the peeling 
properties of the curvature tensor up to order t~ 3 along any forward light cone, but peeling properties 
of higher order require stronger decay of initial data. 

Acknowledgments: The authors would like to thank Sergiu Klainerman for valuable suggestions 
and discussions. 

2 Strategy of the proof and outline of the paper 

2.1 Notations and conventions 

Coordinates: 

• {^ a }a=o,..,3 = (t,x) with t = x° and x = (x\, X2, X3), r = \x\ are the standard space-time 
coordinates 

• s = r + t,q = r — t are the null coordinates 
Derivatives: 

• V = (<9i,<92,d3) denotes spatial derivatives 

• d = (dt, V) denotes space-time derivatives 

• <^ denotes spatial angular components of the derivatives 

• 8 = (dt + d r , (/)) denotes derivatives tangent to the light cones t — r =constant 

• d s = \(d r + dt), d q = \(d r — dt) denote the null derivatives 
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Metrics: 

• m = —dt 2 + Y^=i(d xl ) 2 denotes the standard Minkowski metric on ]R 3+1 

• g denotes a Lorentzian metric solution of the Einstein equations 

• Raising and lowering of indices in this paper is always done with respect to the metric m, i.e., 
for an arbitrary n-tensor Ti ai ,a 2 ..,a n we define n Q ^ 2 an = m ai ^Hp :a2: .. jQn . The only exception is 
made for the tensor g a/3 which stands for the inverse of the metric g a p. 

Null frame: 

• L = dt + d r denotes the vector field generating the forward Minkowski light cones t — r =constant 

• L = dt — d r denotes the vector field transversal to the light cones t — r =constant 

• Si, S2 denotes orthonormal vector fields spanning the tangent space of the spheres t =constant, 
r =constant 

• The collection T = {L, S±, S2} denotes the frame vector fields tangent to the light cones 

• The collection U = {L,L, Si, S2} denotes the full null frame. 
Null forms 

• Qa^(d(/),dil)) = d a (j)dpip — d a <f>d a ^, Qo(d(f>,dip) = m^d^dpip are the standard null forms 
Null frame decompositions: 

• For an arbitrary vector field X and frame vector U, Xy = X a U a , where X a = m a pXP . 

• For an arbitrary vector field X = X a d a = X L L + X^L + X Sl Si + X s ' 2 S 2 , where X L = 
-X L /2, X± = -X L /2, X s * = X Si . 

• For an arbitrary pair of vector fields X, Y: X a Y a = -X L Y L /2 - X L Y L /2 + X Sl Y Sl + X S2 Y S2 

• Similar identities hold for higher order tensors 
Minkowski vector fields {Z} = Z: 

• Z = {d a , Q a p = x a df3 - xpda, S = x a d a }. 

• Commutation properties: [D,9 a ] = [□, Cl a p] = 0, [□,£] = 2D, czO ■= [Z,0]. 
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2.2 Strategy of the proof 



Our approach to constructing solutions of the Einstein (Einstein-scalar field) equations is to solve the 
corresponding system of the reduced Einstein equations (|1.3p 



It is well known that any solution (g, tp) of the full system of Einstein equations, written in wave 
coordinates, will satisfy (|2.ip . The wave coordinate condition can be recast as a requirement that 
relative to a coordinate system {x M } M =o,...,3 the tensor verifies 



Conversely, any solution of the system (I2.ip with initial data (g^ u \t=o, dtg^ u \t=o) and (ip\t=o,dtip\t=o) 
compatible with ()2,2p and constraint equations, gives rise to a solution of the full Einstein system. 
In addition, thus constructed tensor g^ u will satisfy the wave coordinate condition (|2,2|) for all times. 
Verification of the above statements is straightforward and can be found in e.g. |Waj . 
Construction of initial data 

The assumptions of the Theorem assert that the initial data (<?0y j &0ij > "0O> ^l) verify the global 
smallness condition (jl.9p . We define the initial data (gfj,u\t=o, dtg^ult^i^lt^i dtip\t=o) for the system 
(|2.1|) as follows: 

(2-3) 9ij\t=o = 90ij, 9oo\t=o = -a 2 , goi\t=o = 0, ip\ t =o = tpo 

(2-4) dt9ij\t=o = —2akoij, d t g o\t=o = 2a 3 g l 3 k 0ij , d t ip\ t=0 = aipi, 

(2.5) d t goe = a 2 g l J dj go ie - - a 2 g l J d £ g 0ij -ad e a. 

The lapse function o 2 = (l-Mx(r)r -1 ). The data constructed above is compatible with the constraint 
equations and satisfies the wave coordinate condition (|2.2p . In fact, it simply corresponds to a choice 
of a local system of coordinates {x^}^=o,...,3 satisfying the wave coordinate condition \D g x^ = at 
x° = t = 0. This procedure is standard and its slightly simpler version can be found e.g. in [Wa], see 



The smallness condition -Ejv(O) + M < e for the initial data (goij, koij, ipo,ipi) implies the smallness 
condition 



(2.1) 



□<? 9^ = F ilv {g){dg, dg) + 2d^ip d„ip, 
□ 5 ^ = 



(2-2) 




also [LT12] . 



£ N (0) +M < e. 



for the energy £tv(0): 



£iv(0)= Yl ll(l + 1/2+7+|/| VW|, =0 || L2 + £ Kl + r^^W'dth 





+ E (IK 1 + 1/2+7+|/| V vVl*=olb + 11(1 + r) 1 /2+T+|/| V'dt^oh^) . 



0<|i"|<iV 
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defined for a tensor h\ u which appears in the decomposition 

g = m + h° + h 1 , where hP au = xCj)x(r) — <W 

t r 

This allows us to reformulate the problem of global stability of Minkowski space as a small data global 
existence problem for the system 

(2.6) u g h\ v = F^idh, dh) + 2d^ d v i> - a g h°, 

□ g V = o 

with = x(r '/t)x{v)^-5 liv . Arguing as in section 4 of |L-R2| we can show that the tensor g^it) = 
m + h^if) + h^vit) obtained by solving (I2.6f) with initial data given in (12, 3j) - (12. 5h defines a solution 
of the Einstein-scalar field equations. Moreover, the wave coordinate condition (|2.2[) is propagated in 
time. 

The above discussion leads to the following result: 

Theorem 2.1. A global in time solution of (j2.6p obeying the energy bound £ N (t) < Ce{l + tf for some 
sufficiently small 5 > gives rise to a future causally geodesically complete solution of the Einstein 
equations (jl.ip converging to the Minkowski space-time. 

The proof of the geodesic completeness can be established by arguments identical to those in |L-R2| 
and will not be repeated here. 
The mass problem 

A first glance at the system (|2.ip suggests that a natural approach to recast the problem of global 
stability of Minkowski space as a small data global existence question for a system of quasilinear wave 
equations is to rewrite (|2.ip as an equation for the tensor h = g — m: 

(2.7) u g v = F^idh, dh) + 28^ u i>, 

The initial data for h possess the asymptotic expansion as r — > oo: 

(2.8) Vl*=o = 7V + Oir- 1 '*), d t h^\ t=0 = 0( r - 2 - a ) 

for some positive a > 0. While the data appears to be "small" it does not have sufficient decay rate 
in r at infinity due to the presence of the term with positive mass M. This could potentially lead 
to a "long range effect" problem, i.e., the decay of the solution at the time-like infinity t — > oo is 
affected by the slow fall-off at the space-like infinity r — > oo. In the work |C-K| this problem was 
resolved (albeit in a different language) by taking advantage of the fact that the long range term is 
spherically symmetric and the ADM mass M is conserved along the Einstein flow. Thus differentiating 
the solution (in [C-KJ this means the Weyl field corresponding to the conformally invariant part of the 
Riemann curvature tensor of metric g) with respect to properly defined (non-Minkowskian) angular 
momentum and time-like vector fields one obtains a new field still approximately satisfying the Einstein 
field equations but with considerably better decay properties at space-like infinity. We however pursue 
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a different approach by taking an educated "guess" about the contribution of the long range term 
MSpv/r to the solution. Thus we set 

r M 

h% = x(-)xO)— <W 

split the tensor h = h 1 + h° and write the equation for the new unknown h . The important cancelation 
occurring in a new inhomogeneous term U g h° is the vanishing of U^- =0 away from r = 0. The cut-off 
function x( r A) ensures that the essential contribution of U g h° comes from the "good" interior region 
±t < r < |f - the support of the derivatives of x(r/i). Finally, observe that in the presence of the term 
containing the mass M the tensor h = h° + h 1 has infinite energy £/v(0). 

The system (|2.6p is a system of quasilinear wave equations. The results of |ClJ, [K2J show that a 
sufficient condition for a system of quasilinear wave equations to have global solutions for all smooth 
sufficiently small data is the null condition. However quasilinear problems where the metric depends 
on the solution, rather than its derivatives, as in the problems arising in elasticity, do not satisfy the 
null condition. Moreover, as shown in |CBll ICB2] even the semilinear terms F fJiU (dh,dh) violate the 
standard null condition. 

Connection with the weak null condition 

Consider a general system of quasilinear wave equations: 

(2.9) Ufa = A i*P Qfi $ K + cubic terms 

|a|<l/8|<2, |/3|>1 

The weak null condition, introduced in [L-Rl] requires that the asymptotic system for = r(j)j 
corresponding to (I2.9j) : 

(2.10) (d^X^-d^^r- 1 E A inm(dt-d r ) n ^J {d t -d r ) m <S> K , 

n<m<2, m>l 

has global solutions for all small data. Here, the tensor 

Ainm(u) ■= ("2)- m " n £ AjXpCTCbP i=H, W ), W GS 2 

|a|=n, |/3|=m 

The standard null condition, which guarantees a small data global existence result for the system 
(|2.9p . is that Aj^ m (uj) = and clearly is included in the weak null condition, since in that case the 
corresponding asymptotic system (I2.10p is represented by a linear equation. Asymptotic systems were 
introduced by Hormander, see [Hl| IH2] , as a tool to find the exact blow-up time of solutions of scalar 
equations violating the standard null condition with quadratic terms independent of cj>, i.e. \a\ > 1 
in (|2.9p . This program was completed by Alinhac, see e.g. [Alj . It was observed in |L1| that the 
asymptotic systems for quasilinear equations of the form Ucj) = s ^ j c a ^cf) d a ds<p in fact have global 
solutions. In other words these equations satisfy the weak null condition. It was therefore conjectured 
that these wave equations should have global solutions for small data. This has been so far only proven 
for the equation 

(2.11) U(f> = (t>A<p, 
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|Llj (radial case), |A3] (general case). Note that the case \a\ = \(3\ = 0, in (|2.9|) is excluded. The 
asymptotic system only predicts the behavior of the solution close to the light cones, and for the case 
□(/> = 4> 2 the blow-up occurs in the interior and much sooner, see |J2[ IL3| 

The asymptotic system (|2.10p is obtained from the system ()2.9p by neglecting derivatives tangential 
to the outgoing Minkowski light cones and cubic terms, that can be expected to decay faster. In 
particular, 

\3(f> = r~ 1 (dt + dr)(dt — dr)(r(j)) + angular derivatives, 
dfj, = — \Cj^(dt — d r ) + tangential derivatives 

Recall that for solutions of linear wave homogeneous wave equation derivatives tangential to the forward 
light cone t = r decay faster: □<^ = 0implies that \d(j)\ <C/t while \dft\ <C/t 2 . 

A simple example of a system satisfying the weak null condition, violating the standard null con- 
dition and yet possessing global solutions is 

= 03 • 0Vl + (%) 2 , 

□02 = o, n<f> 3 = o. 

Another, far less trivial example is provided by the equation (|2,11|) , 

The asymptotic system for Einstein's equations can be modeled by that of (|2.12p . To see this we 
introduce a null-frame {L, L, S%, S2} decomposition of Einstein's equation. With h = g — m we have 

(2.13) n g h^ = F flu {h){dh,dh) + 2d fJ ,<iijd v dip, Q s v = o 

where 



F, 



(h)(dh,dh) = \d^d v h p - \d^d v h aP + Q^dh.dh) + Gp,(h)(dh,dh), 



fJ,U w V 1 1 A- p--U-v-p 2 

with Q^ u - linear combinations of the standard null-forms and G fiu (h)(dh,dh) contains only cubic 
terms. 

The asymptotic system for the reduced Einstein equations has the following form: 

(2.14) (d t + d r )8 q D^ = H LL d 2 D^ - {P(d q D, d q D) + 2d q <S> d q $) , d q = d t - 8 r 

(2.15) (d t + d r )d q $ = H LL d 2 q <z>, 

ZV~ r V» H a ^ = g^-m a ^ H LL = H a/3 L a L^ L a d a = d t + 8 r 

Here 

(2.16) P(D,E) = \D«E/ - l -D^E ap . 

On the other hand the asymptotic form of the wave coordinate condition (|2.2p is 

(2.17) d q D LT ~0, T €T = {L,S 1 ,S 2 } 

Observe that (|2.17p combined with the initial asymptotic expansion (|2.8p of the metric g suggests that 
asymptotically, as t — > 00 and \q\ = \t — r\ < C, 

M 

(2.18) h La ~ H La ~ — < 
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Decomposing the system with the help of the null frame {L, S\, £2} = T and L : L a d a = d r — dt and 
using that L^u^ = A^u^ = 0, we obtain that 

(2.19) (d t + d r )d q D L L = H LL d 2 q DLL - 2r~ l P{d q D, 8 q D) - r -%$ d q $ 

(2.20) (d t + d r )d q D TU = H LL d 2 g D TU , TeT={L,S u S 2 }, U eU = {L,L,Si,S 2 }, 

(2.21) (dt + d r )d q <$> = H LL d$$ 

In view of (12.180 the equations (I2.20p allows us to find the asymptotic behavior of the components 
dqDTU ~ const and consequently hru ~ t • The asymptotic behavior <3? ~ can be similarly 
determined from (12.210 . On the other hand, 

P(d q D, d q D) = ^d q D LL d q DLL + d q D TU ■ 3 q D TU 

It is the the absence of the quadratic term (dgDiL) 2 and the boundedness of the components dgDxu 
that allows us to solve the equation (I2.19P for Dll although the suggested asymptotic behavior is 
different H ll ~ Hll ~ hit. 

It turns out that the asymptotic system indeed correctly predicts the asymptotic behavior of the 
tensor g^ u = m^ u + h^ u - solution of the Einstein equations in wave coordinates 

Before proceeding to explain the strategy of the proof of our result we review the ingredients of a 
proof of a typical small data global existence for quasilinear wave equations in dimensions n > 4 or 
equations satisfying the standard null condition. For simplicity we consider a semilinear equation 

(2.22) U<j) = N(d(j),d^) 

with a quadratic nonlinearity N. We first note that almost every small data global existence result 
was established under the assumption of compactly supported (or rapidly decaying) data, which by 
finite speed of propagation ensures that a solution is supported in the interior of a light cone r = t + C 
for some sufficiently large constant C. The proof is based on generalized energy estimates 

(2.23) E N (t) < E N (0) + Yl f W zlF i,r)\\^El{r)dr 

Z£Z,\I\<N 

for solutions of an inhomogeneous wave equation Ucft = F. The generalized energy 

E N (t)= £ wdz'tmb' 

zez,\i\<N 

with vector fields Z = {d a , VL a p = —x a dp + xpd a , S = x a d a }, of a solution of ()2.22|) can be shown to 
satisfy the inequality 

(2.24) E N (t)<exp[C sup I \\dZ J 0(t)\\ l ~}E n {O) 

Z£Z,\I\<N/2+2 JO 

The other crucial component is the Klainerman-Sobolev inequality which asserts that for an arbitrary 
smooth function ip 

n—l 1 

(2.25) \difj(t,x)\ <C(l + t + r)-—(l + \q\)-2 E A (t), r = \x\, q = r - t. 
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Combining (|2.25|) with the energy inequality leads to the proof of the small data global existence result 
for a generic semilinear wave equation in dimension n > 4. It also elucidates the difficulty of proving 
such a result in dimension n = 3. In fact, as was shown in [J1J . the result can be false in dimension 
n = 3, e.g., the equation □</> = (j) 2 admits small data solutions with finite time of existence. It is 
interesting to note that the corresponding asymptotic system for such equation leads to a Ricatti type 
ODE. 

In the case when the semilinear terms N(d(p, deft) obey the standard null condition, i.e., it is a linear 
combination of the quadratic null forms Q a p, Qo, it is possible to refine the energy inequality (I2.24D so 
that a combination of the energy and Klainerman-Sobolev inequalities still yields a small data global 
existence result in dimension n = 3. This can be traced to the following pointwise estimate on a null 
form Q: 

\Q{d<t>, 8<f>)\ <C(l + t + r)- 1 £ \8<f>\ \Z<j>\ 

In the absence of the standard null condition in dimension n = 3 the inequalities (|2.24J) -( j2.25p are just 
barely insufficient for the desired result. An illuminating example is provided by a semilinear version 
of the system fl2~T2|) 

D<h. = (dfo) 2 , Ufa = 

A combination of (|2.24|) - (|2.25p applied to the vector <f) = (<pi, ^2) would incorrectly suggest a possible 
finite time blow-up. This analysis however reflects the following phenomena: small data solutions 
of the system above have a polynomially growing energy Ejy(T) ~ t s and the asymptotic behavior 
d(j)(t) ~ t~ 1 lnt. Of course the small data global existence result in this case by applying (]2.23j) - 
(|2.25p separately to each of the components of (j> = (<fii,<f>2)- Recall however that the system (|2.12|) 
models the Einstein equations only after the latter is decomposed relative to its null frame components 

hjjL, h,L,L Such decompositions do not commute with the wave equation and thus prevent one from 

deriving separate energy estimates for each of the null components of h. In this paper we are able 
to solve this problem by adding another ingredient: an "independent" decay estimate. The discussion 
below will be focused on the tensor h 1 = g — m — h° obtained from the original metric g by subtracting 
its " Schwarzschild part". We reluctantly allow for the fact that the energy Ej\f(t) of the tensor h 1 
might be growing with the rate of e 2 t s as t — > 00 for some small constant 5 > dependent on the 
smallness of the initial data. The Klainerman-Sobolev inequality (12.25H would then imply that 

(2.26) sup \dZ I h 1 (t,x)\<Ce{l + t + r)- 1+5 (l + \q\y^ 

ZeZ,\I\<N-3 

In order to close the argument, i.e., verify that the energy E^it) indeed grows at the rate of at most 
e 2 t s we must, according to (j2.24|) . upgrade the decay estimates (12.261) to the decay rate of et^ 1 . To 
do that we invoke the asymptotic system (|2.14p - (|2.15p . It turns out that merely assuming that the 
energy E]\[(t) < e(l + t) s and consequently ()2.26|) allows us to show that the asymptotic system 
provides an effective approximation of the full nonlinear system, i.e., the discarded terms containing 
tangential derivatives do not influence the asymptotic behavior of the field h . The crucial property of 
the asymptotic system f)2.14|) - (j2. 15[) is that, as opposed to the full nonlinear equation, it does commute 
with the null frame decomposition of h thus leading to the system (|2.19p - (j2.2ip . Once the validity 
of the asymptotic system is established the asymptotic behavior of its solutions provides a new decay 
estimate and, as we described above, shows that for T = {L, Si, S2}, U = {L, L, Si, S%} 

\dh\ ru < CeT 1 , \dh\cc < Cet^hxt, \d^\ < Cet~ l 
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The decay rate of et~ 1 lnt can be potentially disastrous since (|2.24|) would imply that the energy 
En(£) could grow at the super-polynomial rate of e 2 exp[Celn 2 1]. However the remarkable structure 
of the Einstein equations comes to the rescue. The analysis of the semilinear terms F„ v (dh, dh) in the 
equation 

U g h\ v = F^dh, dh) - U g h% + 2d^ d v if> 

shows that the sharp t^ 1 decay is only required for the dhj-y components in order for the energy 
estimate -EW(i) < Ce 2 t s to hold. A glimpse of this structure has already appeared in our discussion of 
the weak null condition. Recall that 

F^(dh, dh) = p^h, d u h) + Q^idh, dh) + G^(h)(dh, dh). 

The cubic term G and the quadratic form Q satisfying the standard null condition are consistent with 
the uniform boundedness of the energy E]y(t) and thus irrelevant for this discussion. The quadratic 
form P(d^h, d u h) can be decomposed relative to the null frame and obeys the following estimate 

|P(^, 3*^)1 < \dh\ru\dh\uu, 
i.e., the most dangerous term \dh\ 2 ^ c , which would lead to the damaging estimate 



E (t)<exp[C f \dh(r)\cc]E (0), 
Jo 



is absent! 

The picture painted above is clearly somewhat simplified: we have only indicated the argument 
implying that the lowest order energy E (t) < Ce 2 t 5 and we have not even begun to address the effect 
of the quasilinear terms. In what follows we describe the building blocks of our result and explain 
challenges of the quasilinear structure. 
Energy inequality with weights 

(2.27) / m 2 w+ f [ \d<P\ 2 w' <8 / \d<f>\ 2 w + Ce [ /* + 16 / [ \n g <j>\\d t <f>\w, 

Jt, t Jo Jut J~£o Jo JT, t i + t Jo Jut 



where d denotes derivatives tangential to the outgoing Minkowski light cones q = r — t. The weight 
function 

'l + {l + \q\) 1+2 ~<, q>0 

i + (i + kir 2/ \ q<o 



w(q) 



with n>0, 7>— 1/2 serves a double purpose: it generates an additional positive space-time integral 
giving an a priori control of the tangential derivatives d(j) (we mention in passing that the use of such 
space-time norms leads to a very simple proof of the small data global existence result for semilinear 
equations □(/> = Q(d(j), d(p) satisfying the standard null condition. To our knowledge this argument 
has not appeared in the literature before), it provides the means to establish additional decay in q via 
the Klainerman-Sobolev inequality. The energy estimate (|2.27p is established under very weak general 
assumptions on the background metric g a P = m a ^ + H a ^, e.g., 



1 , ,s 1 



\dH\ <Ce(l + t)-2(l + |g|)-2(l + g_)-^ q_ = \ mm(q, 0)|, 



14 



which are consistent with our expectations that if the energy E^(t) < e 2 (l + t) s then the tensor h 
and consequently H = —h + 0(h 2 ) decay with a rate of at least t~ 1+s . However, as is the case with 
virtually every estimate in this paper, special stronger conditions are required for the Hll, component 
of the metric: 

\dH\ LL <Ce{l + t)-\l + \q\)-\ 

Once again we note that such decay is consistent with the behavior predicted by the asymptotic 
system for the Einstein equations, yet even better control for this component is provided by the wave 
coordinate condition, which we will turn to momentarily. 

The interior estimate (|2.27p . i.e., with w(q) = for q > 0, in the constant coefficient case basically 
follows by averaging the energy estimates on light cones used e.g. in (SJ. We also note that the 
interior energy estimates with space-time quantities involving special derivatives of a solution were 
also considered and used in the work of Alinhac, see e.g. |A2j . [A3] ) . In |L-R2j we proved the interior 
estimate (|2.27p under natural assumptions on the metric g following in particular from the wave 
coordinate condition. The use of the weights w(q) in the exterior q > in energy estimates for the 
space part J" s \d<ft\ 2 w originates in |K-N2j . 

Motivated by the energy inequality (|2.27p we define the energy Ejf(t) associated with the tensor 

r M 

solution of the reduced Einstein equations r} g h = F^ u — D g h° + d^tpd u tp and a collection Z = 
{d a , Q a /3 = —x a dp + xpd a , S = x a d a } of commuting Minkowski vector fields, as follows: 

E N (t) = Yl (ll^ 1 dZ'hHt, .)\\\ 2 + dZ^(t, -)||! 2 ) 
zez,\i\<N 

The Klainerman-Sobolev inequalities with weights 

v, 1 /* QdZ'hHt, x)\ + \dZ^(t, x)\) < h ^ E " { ^ hl/2 , \I\<N-2 

(l + i + rj^l + lgl) 1 /^ 

is the fundamental tool allowing us to derive first preliminary, sometimes referred to as weak, decay 

estimates. 

Decay estimate 

The following decay estimate 



2.28) w( g )(l+t)|^(*,x)|</V+T)||w5^(r,0IU«dr+V [ \\zvZ 1 ftr, -)\\ L 

Jo m2 Jo 



dr 



is the additional "independent" estimate designed to boost the weak decay estimates derived via the 
Klainerman-Sobolev inequality. The weight function 



w(q) 




q>0 
q < 



with 7' > — 1,// < 1/2 is chosen in harmony with the weight function w(q) of the energy estimates. 
The unweighted version of this estimate was used in [LI] in the constant coefficient case. We generalize 
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this estimate to the variable coefficient operator D g = g a ^d a dp under very weak general assumptions 
on the metric g a ^ = m a P + H a ^ , e.g., 



f\\H{T,-)\\ L ^^<e. 
Jo 1 + T 

However, once again the Hct components are required to obey the stronger condition 

\dH\ CT <Ce{l + ty\l + \q\)- 1 

The estimate above is closely connected with the asymptotic equation in the sense that it is obtained 
by treating the angular derivatives as lower order and integrating the equation 

{d t + d r )d q {rcl)) = rU4>+\A^ 

An important virtue of (I2.28P is that when applied to systems = F^u it can be derived for each 

null component of (p. This property is indispensable in view of the fact that the weak null condition 
for the system of reduced Einstein equations becomes transparent only after decomposing the tensor 
h^v relative to a null frame {L,L, A, B}. In particular, (|2.28|) will lead to the estimates 

\dil>\ + \dh\ru < Ce{\ + t + \q\y\ \8h\ < Ce{\ + ty 1 ln(l + t). 

Commutators 

As was mentioned above the combination: "energy estimate - Klainerman-Sobolev inequality- 
decay estimate" , results in sharp decay estimates for the first derivatives of the tensor h 1 and merely 
polynomial growth of the energy Eo(t). However, both the Klainerman-Sobolev and the additional 
decay estimate require control of the higher energy norms En (t) involving the collection Z of Minkowski 
vector fields. To obtain that control one needs to commute vector fields Z £ Z through the equation 
Ogh 1 = F — O g h° and reapply the energy estimate (|2.27p . This is by far the most difficult task in 
dealing with quasilinear equations. The collection Z enjoys good commutation properties with the wave 
operator □ of Minkowski space: [Z, □] = — c^D, where cz = 2 for Z = S and is zero otherwise, however 
its commutator with the wave operator \D g can create undesirable terms. It is exactly for this reason 
that both in the work |C-K| on stability of Minkowski space and |A3| on small data global existence for 
the equation Ucft = (f>A(j) a different collection of geometrically modified vector fields Z was used. The 
new vector fields are adapted to the true characteristic surfaces of the metric g (this originates and is 
especially manifest in the beautiful construction of |C-Kj ) and better suited for commuting through 
□ 9 , this construction however adds another complicated layer to the proof. 

In this work we employ the collection Z of original Minkowski vector fields and argue that the bad 
commutation properties of Z with \D g can be improved if one takes into account the wave coordinate 
condition 

(2.29) djg^y\detg\) =0 



satisfied by the metric g in a coordinate system {x /i } A1= o,...,3. To explain this phenomenon consider the 
commutator [Z, \D g ] with one of the Minkowski vector fields Z G Z. For simplicity assume that Z ^ S 
so that [Z, □] = 0. Then 

[Z, D g ] = [Z, H a/3 d a dp] = (ZH a ? + ntf) d a dp, Hf := H^c^ + H^ C y 
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and the coefficients cj = d a Z@ . It turns out that since Z is either Killing or conformally Killing vector 
field (Z = S) vector field of Minkowski space the coefficients c£ have the property that = 0., which 
implies that \H^r L \ < \H\ct- I n accordance with the usual arguments the worst term generated by 
the commutator [Z, C\ g ] contains two derivatives d q transversal to the light cones q = r — t and any 
modification to the vector fields Z targets to eliminate such a term. In our case this term comes with 
a coefficient (ZH— + H^) and thus can be estimated by 

ZH\r r + \H\ ct , _ „, i , 



1 + M 

Examining the equation 

UgZh 1 = [D g , Z\h l + ZF - ZOgh + Z(dip dip) 

it is not too difficult to see that in order to control (i.e. establish a polynomial bound) the energy 
E\{t) one needs to at least show that 

\ZH\cc + \H\cr < Ce(l + ty l (l + \q\) 

To explain the difficulty of such estimate we note that in the region \q\ < C this estimate is saturated by 
the Schwarzschild part x( r /t)x{ r )M5^ u /r of H. While the desired estimate for Hct can be obtained 
from the sharp decay estimates on the first derivatives of h 1 , the estimate for ZHll is a much more 
subtle issue. An attempt to return to the asymptotic system for the Einstein equations or alternatively 
the decay estimate (|2.28p and show that ZHll still satisfies the sharp decay estimate requires once 
again a commutator argument, this time in the context of the decay estimates, and fails. In fact, we 
can only show that the higher Z derivatives of h 1 decay at the rate of Ce(l + t)~ l+Ce . It is at this 
point that we recall that the tensor g a P = m a ^ + H a ^ verifies the wave coordinate condition (j2.29|) 
which after simple manipulations imply that 

(2.30) dqHcr = dH + 0(H ■ OH) 

Here 5 stands for a tangential derivative, which in particular means that \dH\ < (1 + t)~ l \ZH\. Thus 
if we expect H and its higher Z derivatives to decay with the rate i _1+<5 consistent with the t s energy 
growth and the Klainerman-Sobolev inequality we can conclude that d q HcT w iU decay with the rate 
of at least t~ 2+s . Integrating with respect to the q = r — t variable from the initial data at t = we can 
get the desired rate of decay for Hct- We should note that in order to exploit the additional decay in 
t we need to keep track and use the decay in \q\. This, in particular directly applies to the integration 
described above. The necessity of utilizing the decay in q explains our desire to work with weighted 
energy and decay estimates. Amazingly (|2.30p is preserved after commuting with a Z vector field, but 
only for the Hll component, i.e., 

d g ZHcc = dZH + ZO(H ■ dH). 

Moreover, further commutations would destroy this precise structure. Luckily, the control of higher 
energies Ek{t) still only requires the estimate 

\ZH\cc + \H\ct < Ce(l + t)- x (l + \q\), 
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since the principal term in the commutator [Z^,D 9 ] is still (ZH— + Htt~ )d^Z I . We should also 
mention that from the point of view of the energy estimates the commutator [Z 1 , Cy/i 1 contains 
another dangerous term: 

Since most of the Z derivatives fall on the fist factor one is forced to take it in the energy norm and use 
the decay estimates for d^h 1 . Note the two essential problems: the term Z 1 Hcc does not contain a 
derivative d, which is required in order to identify with the energy norm, the decay of d^h 1 falls short 
of the needed decay rate of t . We resolve these problems as follows. First, a Hardy type inequality 
allows to convert a weighted I? norm of Z 1 Hcc m t° a weighted I? norm of its derivatives. We then 
use the wave coordinate condition, which even after commutation of Z 1 preserves some of its strength: 

|9^ J |cr< Yl \9Z J H\+ \8Z j H\ + Z i O{H -dH) 

\A<\i\ \A<\i\ 

The crucial fact here is that the principal term dZ 1 H contains tangential derivatives 8 and thus can 
be compared to the positive space-time integral on the left hand-side of the energy inequality (|2.27p . 
which in this case will be 

f [ {dZWfv/, 
Jo Jt, t 

which means that we no longer have to be concerned about the decay of d^h in t but rather trade it 
for the decay in q, also needed in the Hardy type inequality, which we always have in abundance. 

Scalar field 

The scalar field ip satisfies the wave equation O g ip = and contributes the quadratic term of the 
form d^ip d^ijj to the right hand side of the equation for the tensor h^ u . Its asymptotic behavior is 
determined by the decay rate \dip\ < Ce(l + i+ \q\)~ ■ All of these properties indicate that tp is very 
similar to the "good" components hju of the tensor h. The only substantial difference, which in fact 
works to the advantage of tp, is that unlike hfu function ip satisfies its own wave equation and thus 
admits its own independent energy estimates. This discussion indicates that in the problem of stability 
of Minkowski space in harmonic gauge, coupling the scalar field to the Einstein equations does not lead 
to fundamental changes in the structure of the equations and requires only superficial modification of 
our analysis in the vacuum problem. As a consequence, starting from Section 9, the proof will be given 
only for the vacuum case with tp = 0. 

We conclude this section by giving the plan of the paper. For the sake of simplicity and clarity of 
the exposition we will only deal with the case of a vanishing scalar field. The generalization to the case 
of non-vanishing scalar field is at each point immediate. Theorem 11.11 for the case of vanishing scalar 
field is stated in section [9j where the actual proof of the theorem starts. 

In Section 3 we write the Einstein equations as a system of quasilinear wave equations for the 
tensor h^ u = — m M „ relative to a system of wave coordinates {x At } At= o,...,3- 

In Section 4 we define the Minkowski null frame {L, L, Si, S2} and describe null frame decomposi- 
tions. We estimate relevant tensorial quantities, including the H\ g and the quadratic form P appearing 
on the right hand side of the reduced Einstein equations, in terms of their null decompositions. 
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Section 5 introduces the collection of Minkowski vector fields Z and records relations between the 
standard derivatives d and vector fields Z. It contains an important proposition giving the estimate 
for the commutator between D g and the powers Z l . More details are contained in the Appendix A. 

In Section 6 we derive our basic weighted energy estimate. 

Section 7 deals with the decay estimates for solutions of an inhomogeneous wave equation O g 4> = 
F. The key result is contained in Corollary 17.21 

In Section 8 we discuss the wave coordinate condition. In particular, we state the estimates for 
the Hcc an d HcT components of the tensor H a/3 = g a ^ — m a ^ and its Z derivatives. Some of the 
details are provided in Appendix D. 

Sections 9-11 contain the statement and the proof of the small data global existence result for 
the system of reduced Einstein equations. 

In section 9 we state our main result and set up our inductive argument. We assume that the 
energy £/v(^) of the tensor fa* = — h® obeys the estimate i?7v(i) < Ce 2 (l + t) s and on its basis 
derive the weak decay estimates for h and h 1 and estimate the inhomogeneous terms F^ v , F® = O g h^ u 
appearing on the right hand side of the equation for h . The weak decay estimates are obtained by 
means of the Klainerman-Sobolev inequality proved in Appendix C. 

In section 10 we use the decay estimates derived in section 7 to upgrade the weak decay estimate 
of section 9. 

Section 11 uses the energy estimate of section 6 to verify the inductive assumption on the energy 
£jv(i). The results of section 11.3 heavily rely on Hardy type inequalities established in Appendix B. 



3 The Einstein equations in wave coordinates 

For a Lorentzian metric g and a system of coordinates {X 1 *} u=o,..,3 we denote 

(3.1) T x v = -g XS '(dpgsv + dvd&v ~ d s9^u), 

the Christoffel symbols of g with respect to {x^}. We recall that 

(3-2) R^vS = ds^^v — d u T 'J's + r^T^ — T p x v T ? 5 

is the Riemann curvature tensor of g and = RJ* va is the Ricci tensor. 

We assume that the metric g together with a field ip is a solution of the Einstein-scalar field 
equations 

(3.3) R^y = d^ijj d u ifj. 

A system of coordinates {x^} is called the wave coordinates iff the Cristoffel symbols T verify the 
condition: 

(3.4) ^1^ = 0, VA = 0,..,3 

Each of the following three equations is equivalent to (|3.4p , which we will refer to as the wave coordinate 
condition, 

(3.5) WW) = 0, g aP d a g^ = \^d^ d a g™ = l -g a pg^d,g^ 
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The covariant wave operator U g in wave coordinates, in view of (|3.5|) . coincides with the reduced wave 
operator g = g a ^d^, 

(3.6) U g = D g = -^=d a g a ? y/\f\d p 

V\9\ 

Define a 2-tensor h according to the decomposition 

where m is the standard Minkowski metric on M 3+1 . Let m^ u be the inverse of m^„. Then for small h 
H nu ._ g »u _ m nv = _ h »v + o^fa 2 ), where hS™ = mP" ' m uv 'h^ 

and 0^ v {h 2 ) vanishes to second order at h = 0. Recall that according to our conventions the indices 
of tensors h^ y ,H^ v are raised/lowered with respect to the metric m. The following proposition was 
established in the vacuum case in |L-R 1] and amounts to a rather tedious calculation. 

Proposition 3.1. Let (g,ip) be a solution of the Einstein-scalar field equations then relative to the 
wave coordinates the tensor h and ip solve a system of wave equations 

(3.7) n g h„ u = F^(h)(dh,dh) + 2d^d u ^, 5 9 v = o 

The inhomogeneous term has the following structure: 

(3.8) F^ u (h){dh, dh) = p(d^h, d u h) + Q^(dh, dh) + G^(h)(dh, dh) 

(3.9) P(d^h, d u h) = ^m aa ' <9 M /w mP^' d v hpp - -m^ffl^S^^ d u h a >p>, 

quadratic form Q^v(U, V) for each (fi, v) is a linear combination of the standard null forms m a ^d a UdgV 
and daUdpV — dpUd a V and G fJ , l ,(h)(dh,dh) is a quadratic form in dh with coefficients smoothly 
dependent on h and vanishing for h = 0: G liu (0)(dh, dh) = 0, i,e., G is a cubic term. 



Observe that the quadratic terms in (|3.9p do not satisfy the classical null condition. However the 
trace ixh = m^h^ satisfies a nonlinear wave equation with semilinear terms obeying the the null 
condition: 

g^dadpm^h^ = Q(dh, dh) + G(h)(dh, dh) + Q(di>, d^). 
4 The null-frame 

At each point (t, x) we introduce a pair of null vectors (L,L) 

L° = l, L i = x i /\x\, 1 = 1,2,3, and L° = 1, V = —x i /\x\, « = 1, 2,3. 

Let S\ and S2 be two othonormal smooth tangent vectors fields to the sphere S 2 . Si, £2 are then 
orthogonal to the normal u = x/\x\ to S 2 and (L,L, Si, £2) form a nullframe. 
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Remark 4.1. The null frame described above is denned only locally. Replacing orthonormal vector 
fields Si,S2 with the projections 

X 

(4.1) ^= di - UiCj J dj, w = -— 

of the standard coordinate vector fields would define a global frame. Set do = L a d a and di = for 
* = 1,2,3. Then {do,..., 83} span the tangent space of the forward light cone. 

We raise and lower indices with respect to the Minkowski metric; X a = m a pX^ . For a vector 
field/one form X we define its components relative to a null frame according to 

X a = X L L a + X L L a + X A A a , X a = X L L a + X L L a + X A A a . 

Here and in what follows A,B,C denotes any of the vectors S\,S2, and we used the summation 
convention; 

X A A a = X Sl Sf + X S2 S% 
The components can be calculated using the following formulas 

\rL ^ \r V~L ^ \r vA V" 

= — 7, A L, ^-=—^L, = 

where 

Xy = X a Y a . 

Similarly for a two form tt and two vector fields X and Y we define 

itxy = Tr a pX a YP. 

The Minkowski metric m has the following form relative to a null frame 

m LL = rriLL = m LA = m LA = 0, m L L = m LL = -2, m AB = 5 A B, 

i.e. m a pX a Y^ = -2{X L Y^ + X^Y L ) + 8 AB X A Y B , where 5 AB X A Y B = X Sl Y Sl +X S2 Y S K The 
inverse of the metric has the form 

LL LL LA LA n LL LL /o AB cAB 

m = m — = m = m- =0, m - = m- = — 1/2, m =0 
We also define the tangential trace 

(4.2) tik = 5 AB k AB , 
and record the identity 

^d a = TT^L + TT^L + TT A ^ A = ~7r/L - ^7r/L + 7r/A 

We introduce the following notation. Let T = {L, Si, S2}, = {L, -L, £1, S2}, £ = {£} and 
S = {S\, 52}. For any two of these families V and W and an arbitrary two-tensor p we define 

(4.3) \p\ vw = £ \Pth vPw7 \> 

VeV,WeW, 

(4.4) \d P \vw= £ K^)^^^!, 

E/ew,vev,wew, 

(4.5) |a P |vw= £ |(^/3 7 )T Q ^^| 

Ter,vev,wew, 
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Note that the contractions with the frame is outside the differentiation so they are not differentiated. 

Let Q be one of the null quadratic form, i.e. Q a /3(d(f>, dip) = d a <f)dpip — dpfidatp if a ^ (3 and 
Qo(d(p,dip) = m a ^d a (f)dpip. Motivated by (|3.9|) we define the quadratic form P 

(4.6) P(tt, 9) = l -m aa 'mW\ aP e a ,p, - \m^m a ^\ aP 6 a , p , 

The proof of the following result is a simple exercise and we leave it to the reader. 
Lemma 4.2. For an arbitrary 2-tensor tt and a scalar function (j) 

(4.7) \n^B a ^B^\ < (|vr| ££ |^| 2 + |tt| 1^11301), 

(4.8) {Lan^d^l < (\ir\ccW\ + M \8<f>\), 

(4.9) \{B a ^)3^\ < (\dir\cc + \dn\) \dcf>\ + \Bn\ \B<f>\ 

(4.10) K P d a d^\ < {Hcdd 2 ^ + |tt| \dd<l>\) 

For the quadratic form P defined in (|4.6p and a null form Q 

(4.11) |P(vr, 0)| < |tt \tu\0\tu + k \cc\0\ + |tt | \6\cc, 

(4.12) ^{d^dn^mm + mm 

For a function </>, \B(p\ = J2a=o \^a<t>\ > where <9 Q the tangential derivatives defined after (j4. 1[) . is equiv- 
alent to Y1t&t \T a d a <t>\- However \B 2 (f>\ = Y2a/3=o \BaBp(j)\ is not equivalent to T€ -j-\T a S^B a Bp4>\, 
the difference is of order \d<j)\/r which is a main term. Furthermore \B 2 4>\ need not even be equivalent 
to J2sTeT \T a da{SPdp(j))\, the difference is however bounded by \B<j)\/r which is a lower order term. 

Lemma 4.3. For an arbitrary symmetric 2-tensor 7r a/3 and a function (f> 

(4.13) \n al3 d a dp<f) - v LL B 2 q <\> - 2-K LL d s B q <t> - r~ l tr"vr d q <j)\ < \x\jcr\dd<f>\ + |vr| (\B 2 <f>\ + r^B^). 
Proof. Using that B q , B s derivatives commute with the frame {L,L, Si, S2} we obtain that 

^B a B^ - itllB 2 ^ - 2ir Lk d s d q <j) = ttllB 2 ^ + 2irLAd a (A%<f>) - ir LA L° AB^ + t^abA'B^ BiB^ 
Furthemore, 

TTABA'B^BiBj 4> = ttabA'B^Bj <j>) - TVAB^B^Bj <j> 
Decomposing with respect to the null frame we obtain 

(A%&)dj = (A i B i B) L L + (A i B i B)-L + (A i diB) a C 
Note that \A L BiB^\ < Cr~ 1 , since B^ are smooth functions of to = x/\x\ £ S 2 . Now, since B^tOj = 0; 

2(A i B i B)^ = -(A i 3 i B j )u j = A^^Bmjj = AVB*-^- - w f coj) = -d^A^i = -5 AB 

is precisely the null second fundamental form of the outgoing null cone t — r =const. The lemma now 
follows since as pointed out before the lemma 

£ \T a B a {S?B^)\ <C{\8 2 ct>\+r- l \dcl>\). 
S, Te T 

□ 
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Corollary 4.4. Let 4> be a solution of the reduced wave equation O g 4> = F with a metric g such that 

H aH _ g ap _ m ai3 satis fi es ^ e con dition that \H LL \ < \. Then 

(4.14) 

4 ^~^^~ tr fg^ L ' )^ (r ^ )+ 2^l| ~ r \ A ^\ + \ H \LT (r \Bdcf>\ + \dcf>\) + \H\ (r \8 2 4>\ + \d<t>\+r- 1 



where A w = A = 5 lJ didj. 
Proof. Define the new metric 

xafi _ 9 
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The equation g a ^d a dp(f) = F then takes the form 

F 
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-2g L L 

which also can be written as 

U<p + (g a P - m a/3 )d, 
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Let ir a P be the tensor ir a P = (g a/3 — m a ^) Observe that 

Tr a P = (_2/^)- 1 (g a ? + 2m afi g LL ) = (-2g LL y x (H a ? + m a \2g LL + 1)) 

= (-2g L ±)- x {H a P + 2m al3 H L ^) 

Thus, 

(4.15) TT L L = 0, TT Lr = (^g^)" 1 H LT , t?VT = (^Z^)" 1 (t? # + H L l) 

Moreover, | vr | < \H\, since g L — = H L — — ^ and by the assumptions of the Corollary \H L —\ < \. 

Now using ()4.13p of Lemma 14.31 with the condition that till = 0, together with the decomposition 

U(j) = -d 2 t (t> + A(f>=^(d t + d T ){d r - d t )rcf> + = ~d s d q r<() + A w 0. 

we find that the identity □</> + TT a ^d a d^<j) = (— 2<? i -)~ 1 i ? leads to the inequality 

\4d s d q r<p + rir L Ld*<f> + tiir d q <j) + (2g LL )- l rF\ < r|A^| + r\n\ LT \Bd^\ + \n\ (r|<9 2 0| + \B(f>\) 
Finally, identity ()4.15p and a crude estimate \tt\ < \H\ yield the desired result. □ 

5 Vector fields 

The family of vector fields 

Z = {d a , tt aj3 = -x a dp + xpda, S = td t + rd r } 

plays a special role in the study of the wave equation in Minkowski space-time. We denote the above 
vector fields by Z c with an 1 1-dimensional integer index t = (0, ...1, .., 0). Let I = (ix, tfc), where 
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= 1, be a multi- index of length \I\ = k and let Z 1 = Z Ll ■ ■ ■ Z Lk denote a product of k vector fields 
from the family Z. By a sum I\ + I2 = I we mean a sum over all possible order preserving partitions 
of the multi-index / into two multi-indices I\ and I2, i.e., if / = {t\, l^), then I\ = (^ x , tj n ) and 
h = (k n+ i, ■■■, L i k ), where i\, ...,ik is any reordering of the integers 1, k such that i\ < ... < i n and 
i n+ i < ... < ik- With this convention Leibnitz rule becomes Z I (fg) = ^2j 1+ j 2= j(Z 11 f)(Z l2 g). 

We recall that the family Z possesses special commutation properties Z: for any vector field Z G Z 

[z,n] = -c z o, 

where the constant cz is only different from zero in the case of the scaling vector field C5 = 2. We also 
record the following expressions for the coordinate vector fields: 

(5.1) d t =~ 



(5.2) d r = Lo'di 



t 2_ r 2 ' 

tuj l £loi — rS 



t 2 



_ -x j Qjj + ttlpj - XjS _ _ XjS XjX j fl j fioi 
{ ' 1 ~ t 2 -r 2 ~ t 2 - r 2 + t(t 2 - r 2 ) + t 

In particular, 

(5.4) Bo = d s = \{d t + d r ) = ^^f, * = = ^i= -^ + "« . 

Recall that 8 denotes the tangential derivatives, i.e., Span{<9o, 81, 82, 83} = Span{(9 s , Si, S2}. 
Lemma 5.1. For any function (ft and a symmetric 2-tensor ir 

(5.5) (1 + t + \q\)\8ct>\ + (1 + \q\)\d<j>\ <C^2 \ ZI <PI 

1*1=1 

(5.6) \B 2 c t) \ + r~ l \8(P\<-Y. nSri i ' where l^l 2= E l^0| 2 , 

|/|<2 11 a,/3=0,l,2,3 

(5.7) i^a^i < ^(r^g + ^qf|) E 1^1 

Proof. First we note that if r + i < 1 then ()5.5|) holds since the standard derivatives <9 Q are included 
in the sum on the right. The inequality for \8<p\ in (15. 5|) follows directly from (|5.4j) . The inequality for 
1 5/ 1 in (|5.5|) follows from (|5.ip and the first identity in (|5.3p . 

The proof of (15. 6p follows immediately from ()5.4p and the inequality \diUj\ < Cr^ 1 . The inequality 
(|5.7p follows from Lemma 14.31 (|5.5p and the commutator identity [Z,di] = cfd a , which holds with 
constants cf . □ 



Next we state the result following from the above lemma and Corollary 14.41 
Lemma 5.2. Let cf) be a solution of the equation O g 4> = F with g a @ = m a ^ + H a/3 . Then 



(5 . 8) (44 _ _ *^±^) W ) + - J- < (1 + '-^l + m yi J2 \z< 



2g L L " 2g L ± r ) ^ T ' 2g 



r\H\ 



1 + 



|I|<2 
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Proof. By Corollary 



Ad s ~ ?zfTd q - ^—ff— — )d 9 (r<f>) + 



2g L ^ q 2g L ^r J qK YJ ^ 2g^ 

< r|A w 0| + \H\ LT (r\5d<f>\ + \d</>\) + \H\ (r \8 2 (f>\ + \8<f>\ + r _1 M) 

where A w = 5 l:) didj. Here all the the derivatives can be reexpressed in terms of the vector fields Z and 
d q using l5TTj yielding the expression ()5.8j) . Note that 



\dd<f>\ 



< E|/|=il^l E|/|<i I^VI E|j|<2 



l + ~ l + t + ~ (l + |g|)(l + * + |g|) 

□ 



The last result of this section records an important statement concerning commutation between 
the reduced wave operator d fl and the family of Minkowski vector fields Z. As we already explained 
in Section 2 our small data data global existence result for the system of reduced Einstein equations 
O g h l = F — Dgh° is based on controlling the energy and pointwise norms of the quantities Z 1 ^ with 
vector fields Z G Z - the family of Minkowski vector fields. The above control is achieved via the 
energy and decay estimates for solutions of the inhomogeneous wave equation 

UgZ 1 ^ = Z^gh 1 ~ UgZ 1 ^ + Z J F ~ ^Ugh^, Z=Z + C Z 

and therefore requires good estimates on the commutator Z I D g — DgZ 1 . To underscore the importance 
of this commutator we recall that the only two known examples of quasilinear hyperbolic systems, with 
the metric dependent on the solution rather than its derivatives, possessing small data global solutions 
required the use of modified vector fields primarily due to the lack of good estimates for the commutator, 
[A3], |C-K] . The proof of the proposition below together with other commutator related statements is 
contained in Appendix A. 

Proposition 5.3. Let O g = □ + H a/3 d a ds. Then for any vector field Z 6 Z we have with Z = Z + cz 
(5.9) \UgZ^ - &U g 4>\ < 1 Yl E \ ZJR \ \ QZK ^\ 

\K\<\I\, |J|+(|A-|-1) + <|/| 



1 



£ ( £ \Z J H\ LL + £ \Z*H\ LT + £ \Z J "H\)\dZ 

K\<\I\ |J|+(|JC|-1)+<|J| |J'|+(|«-|-1)+<|J|-1 |J"| + (|K|-l) + <|/|-2 



1 + 

where (\K\ - 1)+ = \K\ - 1 if \K\ > 1 and (\K\ - 1)+ = if \K\ = 0. 



K . 



6 Energy estimates in curved space-time 

In this section we establish basic weighted energy identities and estimates for solutions of the inhomo- 
geneous wave equation 

(6.1) U g <t> = F 
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The weights under consideration will depend only on the distance q = r — t to the light cone t = r and 
will be defined by a function w 

W = w(q) > 0, w'{q) > 

The weights will serve a two-fold purpose: 1) to provide additional decay of (j) in the exterior region 
q > 0, which will follow from the energy estimates via a Klainerman-Sobolev type inequality 2) to 
establish an additional a priori bound on a space-time integral involving tangential derivatives d<j) of 
the solution. 

Lemma 6.1. Let (j) be a solution of the equation (|6.ip decaying sufficiently fast as \x\ —> oo. Assume 
that the background metric g is such that the tensor H a @ = g a/3 — m a ^ satisfies \H\ < \. Then with 
oj = x/\x\ 

(6.2) 

(|d t 0| 2 + \V<P\ 2 )w(q)dx + 2 I 2 I \d<P\ 2 w'(q)dxdt < 4 ! (\d t <P\ 2 + \V<f>\ 2 )w(q) dx 

(6.3) + 2 f 2 f \{2d a H aP )dp(j)dt(l) - (d t H a)3 )d a 0d^ + 2Fd t <P\w{q)dxdt 

+ 2 1 2 I \H al3 d a (f)d l3 (j) + 2(u; i H l13 - H 0l3 )d l3 (j)dt(j)\w'(q)dxdt 

Proof. Let <pi = di<j), i = 1,2,3, and (f>t = dt4>. Differentiating under the integral sign and integrating 
by parts we get 



(6.4) 
d_ 
dt 



-g 00 <P 2 t +g l ^ i cP j )w(q)dx- J 2dj(g°^ 2 w(q))dx = 2 J wiq^-g^^tt+g' 3 <Pi4>tj-2g Qj 4>t4>t 3 ) dx 
+ j w{q) ( - (dtg 00 ^ 2 + (dt&hh - 2(d^)<p 2 ) - w'iq) ( - g°°$ + g^fa + 2^$) dx 
= 2 f w(q)( - g 00 (t>t<t>tt - g ij Ma - Zg^Mtj) dx 
+ J w(q){ - (% OO )0? + (dtgV)^ - 2{d 3 g^)<f ) 2 - 2(d ig ^)Mj) dx 

- J w'(q){- 5 °V? + g ij ^j + 2u jg ^ 2 + 2^,^0,0/) dx 

Hence, 

(6.5) ^ f ( - + g ij ^j)w(q) dx = - f w(q) (2& U g <j> - (d t g a ^ a ^ + 2(8^)^) dx 



dt 



w'(q)(g a U a ^ + 2(oo i9 ia - g 0a )(f>t(f>a dx 



Furthermore, with (fi r = ui l cpi = d T <p and <f>^ = (pi — <jj%4> t = 9$ 
m a0 (f> a (f>p + 2<M^W° - m 0a )4> a = -<g + + 2&(wVt + <t>t) = [4>t + </v) 2 + S^M* = 
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Since \H\ < 1/2 we also have that 
The lemma follows. □ 



We now consider the following weight function: 
(6.6) id = w(q) - 



1 + (1 + |<?|) 1+27 , when q>0 
1 + (1 + \q\)~ 2 ^ when g<0 



for some > 0, 7 > — 1. We clearly have 

(6.7) w' < 4w(l + I ^ I > 1 < 167~V(1 + g_) 2 '\ 

where q- = \q\ for q < and = otherwise. 



Proposition 6.2. Let <ft be a solution of the wave equation (I6.ip with the metric g such that for 

(1 + Mr'lHlii + \9H\ll + \dH\ < Ce'(l + t)~\ 
(6.8) (l + M)- 1 !^! + \8H\ < Ce'(l + t)-5(i + | (Z |)-5(i + g _)-f 



Then for any < 7 < 1, and < e' < 7/C1, we /iave 

(6.9) / \d(j)\ 2 w+ f I \B<j)\ 2 w'<%! \d(/)\ 2 w + W [' [ f C ^ 9 ^ 2 + \F\\d<l>\) w 

JT,t Jo Jt. t Jt.0 Jo Jz t ^ 1 + * 



Remark 6.3. Observe that by the Gronwall inequality the energy estimate of the above proposition 
implies t e growth of the energy. 

Remark 6.4. We recall again that the interior estimate (|6.9p . i.e., with w(q) = for q > 0, in the 
constant coefficient case basically follows by averaging the energy estimates on light cones used e.g. in 
[?]. The interior energy estimates with space-time quantities involving special derivatives of a solution 
were also considered and used in the work of Alinhac, see e.g. [A2], |A3| ) . In |L-R2| we proved the 
interior version of ()6.9|) . The use of the weights w(q) in energy estimates in the exterior q > for the 
space part J" s \dcj)\ 2 w originates in |K-N2j . 

Proof. The proof of the proposition relies on the energy estimate obtained in Lemma 16.11 We decom- 
pose the terms on the right hand-side of (|6.2p with respect to the null frame, using Lemma 14.21 



KdaH^d^M] + \(d t H a P)d a <f>dpct>\ < (\(dH) LL \ + \8H\) |9(/>| 2 + \8H\ \d(/>\ 



Therefore, using the assumptions (J6.8P on the metric g, we obtain that 



(6.10) \2{d a H a P)d^d t 4> - (d t H a P)d a ct>dp<j>\ < —m\' + 



l + V ri (l + |g|)(l + g_) 2 M 
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Decomposing the remaining terms we infer that 

\H a Pd a <f>d 0( j)\ + \L a H af3 (d t (l))d^\ < \H\ LL m 2 + \H\\8ct>\ \d(j>\ 
Once again, using the assumptions (j6.8p . we have 



(6.11) \2H^L a d p cj>d t cj> + H^d a cl>d^\ < e'^f\d<(>\ 2 + e' ™ 



Thus, with the help of ([67 



w+ I f \dcj>\ 2 w'<! \d<j>\ 2 w + e f* f 0P^-w + \d^\ 2 —)+ f i \F\\d(f>\w 



's t Jo Jt, t ^ Jy, Jo is T + * 7 ' Jo Jy, t 

and the desired estimate follows. □ 

7 Decay estimates in curved space-time 

In this section we derive L°° estimates for the first derivatives of solutions of the equation 

(7.i) B^ = f 

These estimates are complimentary to the global Sobolev inequalities derived in Appendix C and 
will provide a way to improve upon the decay estimates derived from the weighted energy estimates 
via global Sobolev inequalities. Estimates of these type were obtained in |Llj in the case of the wave 
equation in Minkowski space-time, i.e., constant coefficient metric g, and |L-Rlj for variable coefficients. 
Here, however, we introduce a weighted version of the L°° estimates deigned to capture additional decay 
in |g| = |r — t\. 

For 7' > — 1, // < 1/2 define the weight 



(7.2) w = w{q) 



;i + M) 1+7 \ when q>0 
'1 + M) 1/2 ~ M ' when q<0 



Lemma 7.1. Let <f> be a solution of the reduced scalar wave equation (17. ID on a curved background 
with metric g. Assume that the tensor H al3 = g a/3 — m a " obeys the estimates 



(7.3) 



r°° dt i 

\H\<J, J o < \H\lt<e'- 



/ kl + i 



+ t+\x\ 

in the region Dt = {x : t/2 < \x\ < 2t}. Then for a = max(l + 7', 1/2 — //) 



(7.4) (l+t)w(q)\d<f>(t,x)\ <C sup V \\w{q) Z^t, 

0<r<t mi 

+C f (e'a\\w{q) d<f>(r, .)|U«+(1+t)||c7(<z) F(r, •)|U-(D T )+ E (i+^IMs) Z 1 ^, Ollx-p.)) ^ 
J ° |/|<2 
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Proof. Since by Lemma 15. II 

(7.5) (1 + \t - r\)\d<t>\ + (1 + t + r)|90| < C {Z 1 



r = |x| 
71=1 



the inequality (|7,4|) holds when r<t/2 + l/2orr>2t — 1. 

By Lemma E2] with = 1/2(9* + <9 r ) and <9 g = l/2(d r — dt), 

(7.6) |(4d s - ^5,)a 9 (r0)| < (l + + I^Qr-- 1 £ + |#| r" 1 |S,(r0)| + r\F\ 

Multiplying by the weight w(q) and using that w'(q) < Cw(q)/(1 + \q\) along with the assumptions 
(17. 3D we obtain 
(7.7) 



(4d s - ^a g Mg)a s (^)| < (M + a J^\ w{q) \ dq{r(})) \ + ^ ^M^Vl + c{t + i M<7 )|f| 



in the region t/2 + 1/2 < r < 2i — 1. Let (r, x(r)) be the integral curve of the vector field d s + 
H—(2g L —)~ l d q passing through a given point (t,x) contained in the region t/2 + 1/2 < r < 2t — 1. 
Observe that, by the smallness assumption on if, any such curve has to intersect the boundary of the 
set t/2 + 1/2 < r < 2t - 1 at (r,y) such that \y\ = t/2 + 1/2 or \y\ = 2r - 1. 

Then along such a curve the function tp := w(q)d q (r(p) satisfies the following equation: 

(7-8) |^| + / 

where 

h = C^ f= ^-Mq)\d q rcf>\ + C(l + t)^q)\F\ + C £ ^MA. 

+ + |/|<2 + 

Thus using the integrating factor e~ H with H = f h(s) ds and integrating along the integral curve 
(t, x(t)) from any point (t, x) in the set t/2 + 1/2 < r < 2t — 1 to the first point of intersection (to, %o) 
with the boundary of the set i/2 + l/2<r<2£ — lwe obtain 

\i/>(t,x)\ <exp(^ \\h(a,-)\\ L ™do^\i;(t ,x )\ + J exp (jf ||M<x, dr) ||/(r', dr', 

with the L°° norms are taken over the set t/1 + 1/2 < r < 2t — 1. 

For the points (to, a>o) such that |xo| = to/2 + 1/2 or \xq\ = 2to — 1 we have by (|7.5|> that 

\4>(to,x )\ < Crw(q)\d q ^\ + Cw(q)\</>\ <CY, 

|/|<i 

The desired inequality now follows from ()7.3|) . which implies that \\h(a, Olli 00 da < \ and the 
inequality 

(1 + t + r)\d<t>\ < C \^4>\ + C\d q (r$)\. 
|/|<i 

□ 
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We now state similar estimates for a system 

(7.9) n<t>^ = 

While it is trivial to extend the estimates of Lemma 17.11 to each of the components of </> MJ/ our interest 
lies in the estimates for the null components of <f>. Contracting (17. 9h with the vector fields {L, L, A, B} 
is far from straightforward. We instead exploit that the null derivatives d s ,d q commute with any 
of the vector fields of the null frame. We assume that the weight function w{q) is as in (j7.2j) and 
q = max(l + 7', 1/2 - //)• 



Corollary 7.2. Let <j)^ v be a solution of the reduced wave equation system (|7.9p on a curved background 
with a metric g. Assume that H a P = g a/3 — m a ^ satisfies 

£ ' f°°n rru Ml dt / £ ' , rrl ^ £ ' M + 1 

" ' ' < — , \ti\LT < 



(7.10) \H\ <-, J \\H(t,-)\\ L oo {Dt) - 



+ t - 4 1 ' 4 l + t + \x\ 



in the region Dt = {x 6 R 3 ; t/2 < \x\ < 2t}. Then for any U, V 6 {L, L, A, B} and an arbitrary point 
x £ D t : 



(7.11) {l + t+\x\)\w{q)dcf>{t,x)\ uv < sup V HwCffJZVCr.OIU- 

°^*|/|<l 

+ / (£%|Mg)|d0(V)|uv||L~+(l+r)|K^ 



\I\<2 

Proof. By Lemma 15.21 for each component we have the estimate 



(7.12) 1(4* - 1^3, - ^l±ftt) 85 (^) + ^\<( 1 + n^ + m y-i Y>Vl 

and since d s and <9 g commute with contraction with the frame vectors we get 

(7.13) | (49, - - t %i^)7 Ww) 4- < (l + ^ 4- W )r-' E 1^1 

The proof now proceeds as in Lemma 17. II □ 



r\H\ 



8 The wave coordinate condition 

The results of previous sections underscore the special role played by the Ull components of the tensor 
_ ga/3 _ m a/3 j n energy and decay estimates for solutions of the wave equation O g 4> = F. 
In this section we explain how the wave coordinate condition on the tensor g^ u provides additional 
information about H^l. 

Recall that the wave coordinate condition for metric g in a coordinate system {x M } /i= o,...,3 takes 
the form 

(8.1) d^y/lfetd) =0. 
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Expressing g^ u in terms of the tensor H^ u we obtain 

5""vTdetg[ = {mT + H^) (l - *tr# + 0{H 2 )) 

Therefore, 

(8.2) d^R^-^m^tr H + 0^ u (H 2 )j = 0, where 0" V {H 2 ) = 0{\H\ 2 ). 

Recall also the family of, tangent to the outgoing Minkowski light cones, vector fields T = {L, A, B}. 
The divergence of a vector field can be expressed relative to the null frame as follows: 

(8.3) fyF" = L^F" - L^F* + A^aF" 
We can now easily prove 

Lemma 8.1. Assume that \H\ < 1/4. Then 

(8.4) \8H\cr < \dH\ + \H\ \dH\ 
Proof. It follows from (JS2J) and ((OJ) that 

(8.5) \L^d{H^ - -m^ tr H) \ < \dH\ + \H\\dH\ 

Contracting with T G T and using that m-fi = gives the desired result. □ 

We now state a generalization of the above result containing estimates for the quantities Z 1 ' Hlt 
with vector fields Z 6 Z our family of Minkowski Killing and conformally Killing vector fields. The 
result is a rather tedious consequence of commuting vector fields Z through the wave coordinate 
condition (|8.ip and we postpone the details of the proof until Appendix D. 

Proposition 8.2. Let g be a Lorentzian metric satisfying the wave coordinate condition (|8.1j) relative 
to a coordinate system {3^1^=0,. ..,3- Let I be a multi-index and assume that the tensor H^ v = g^ v — m^" 
verifies the condition 

\Z J H\<C, V| J| < \I\/2, VZ^Z. 

Then for some constant C 

(8.6) \dZ l H\ LT < C"( \dZ J H\ + \8Z J H\ + ]T \Z h H\\8Z h H\^j 

\J\<W \J\<W-l |^|+IK2]<|J| 

(8.7) |<9Z 7 #| LL <C"( \9Z J H\+ ^ \8Z J H\ + ^ \Z h H\\8Z h H\\. 

\J\<\I\ l^l<|/|-2 |/i| + |/ 2 |<|/|,m>2 

Similar estimates hold for the tensor = g^ v — m^. 
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9 Statement of the Main Theorem and beginning of the proof 

We consider the initial data (g(tv\t=o,dt9nu\t=o) f° r the Einstein-scalar field equations R^ u = d^ipd^, 
constructed in (|2.3|) - (|2,5|) . The spatial part go^ of g^ v \t=o together with the second fundamental form 
koij = —l/2dtgij\t=o satisfy the constraint equations 

- l^ol 2 + (trfco) 2 = ift + |V^o| 2 , V^ 0ij - Vitrfco = Vi dupo- 
By construction the initial data (gfiu\t=o,9tg^ u \t=o) satisfies the wave coordinate condition 



The reduced Einstein equations, written relative to a 2-tensor h^ u = g^ u — take the form 

(9.1) Ug V = Fpj, + 28^ d„t/j, D g ^ = 0, 

(9.2) F^(h)(dh,dh)=P{d„h,d„h) + Q^{dh,dh) + G lJlu (h)(dh,dh), 

(9.3) P{d^h,d u h) := ^trhdvtrh - ^h a ^d u h Q p, , 

where Q is a linear combination of standard quadratic null forms and G is a quadratic term in dh 
vanishing in h at h = 0. We assume that the initial data for h(t) are such that at t = 0, g^j, = m^+h^ 
verifies the constraint equations and the wave coordinate condition. 
Define 2-tensor h\ v {£] 

(9.4) h^t) := V(i) - h%,(t), where h%,{t) = x(r/t)xW^ P 

where x( s ) £ C°° is 1 when s > 3/4 and when s < 1/2. It is clear that we can reinterpret (|9.1j) as 
the equation for a new unknown - 2-tensor h 1 : 



(9.5) Dgh 1 ^ = - U g h% + 2d^ d V l/j 

Set 



(9.6) £ N (t) = Yl [W^dZ'hHt, -)\\ L2 + Ww^dZ'^t, OII^J , 

|/|<7V,Ze.Z 



where 



w 



'l + (l + |g|) 1+2 T, g>0 
l + (l + |g|)- 2 ^, g<0 



where q = r — t and 7,/i > 0. We recall that the initial data for ip is given by ip\t=o = "^o an d 
d t ip\t=o = atpi with a 2 = (1 — x(^)Afr _1 ). 
We now state the main result. 

Theorem 9.1. There exist a constant eq > such that ife<Eo and the initial data (h \t=0i dth}\f = Q, i/^o, ipi) 
are smooth, and obey £n(0) + M < e together with the condition 

(9.7) liminf (1^(0,^)1 + |^o(x)|) 0, 

\x\— >oo 
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then the solution h(t) of the reduced Einstein equations (|9.1|) can be extended to a global smooth solution 
satisfying 

(9.8) £ N (t) < C N e(l + t) C£ 
where CV is a constant depending only on N and c is independent e. 

Remark 9.2. Recall that if the initial data for (g^ u ,tp) verifies the constraint equations and the wave 
coordinate condition then there exists a local in time classical solution [h^ v {t), i/)(t)) of (|9.1[) such that 
9uv(t) = m [iv(t) + h^ v {t) obeys the wave coordinate condition for any time t in a maximum interval of 
existence. 

Therefore, in what follows we shall assume that a local in time solution g^ u {t) = m^ v + 
havii) obeys the wave coordinate condition 

(9.9) d^(tW\detg(t)\) =0 
for any < t < T . 

We also note that the maximum time of existence To can be characterized by the blow-up of the 
energy £jv(i) — >• oo as i — >■ T ~~. 

Remark 9.3. As was explained in Section 2 the proof of the result for the Einstein-scalar field problem 
requires only superficial modifications as compared to the vacuum case. 

Therefore, in what follows all the arguments will be provided only for the vacuum 
problem with if) = 0. 

For the proof we let 5 be any fixed number < S < 1/4 and 5 < 7. We define the time T < Tq to 
be the maximal time such that the inequality 

(9.10) £ N (t) < 2C N e(l + t) s 

holds for all < t < T. Note that by the assumptions of the Theorem T > 0. We will show that if 
e > is sufficiently small then this inequality implies the same inequality with 2CV replaced by Cat 
for t < T. Since the quantity is continuous this contradicts the maximality of T and it follows that 
the inequality holds for all T < Tq. Moreover, since the energy £jsr(t) is now finite at t = To we can 
extend the solution beyond To thus contradicting maximality of To and showing that To = +00. 

The first step is to derive the preliminary decay estimates for h 1 under the assumption (I9.10p . The 
estimate (|9.10p can be explicitly written in the form 

(9.11) J2 \\w{q) ll2 dZ I h 1 {t,-)\\ L 2<Ce{l + t) s , 0<£<7 
zez,\i\<N 

with some sufficiently large constant C. The following result is a consequence of weighted global 
Sobolev inequalities proved in Appendix C. 
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Corollary 9.4. Let h 1 verify (|9.11|) and h° be as in (|9.4|) . Then for i = 0,1, 

(9 12) \8Z'hHt x)\ < { C£{1 + 1 + lql) ' 1+5(1 + kir1 " 5 '' q> ° \I\ < N -2 

where 5' = 5, if i = and 6' = j > 8 if i = 1. Furthermore, 

and 

\Ce(l + t+|<,|)- 2+i (l + k|) 1/2 , 9<0 

Proof. We will only prove the estimates for i = 1 since the estimates for i = follow by a direct 
calculation from the form of h . Estimate (j9.12fl follows from the weighted Sobolev inequality of 
Proposition [TTTJ We claim that 

(9.15) lim \Z I h 1 (0,x)\^0, \I\ < N - 2 

\x\—>oo 

In fact, if |/| = this follows from (|97f|) and (j9"T2|) . and if \I\ > 1 this follows from ([9"?T2"j) , since 
\Z(j)\ < C(l + t+ \x\)\d(f)\. Estimate (j9J3|l for t = follows by integrating (j9J2l) from space-like 
infinity, where (|9.15p hold. 

Remark 9.5. The weighted Sobolev inequality of Proposition 114. ll in fact implies the estimate 

i 8 zw (4 ,x)i<{ ce(i+(+i?ir : (i+i?ir :r (i+f ' i ' ,>o , i/i<»-2. 

1 W " - \Ce(l + t+| g |)- 1 (l + kl)^ 1/2 (l + *) i , <f<0 
In particular, 

(9.16) |a^ J ^(0,x)| < Ce(l + Ixl)- 2 " 7 , 

(9.17) IZ^O.z)! < Ce(l + M)" 1 " 7 , 

The estimate (19.130 for r > t follows by integrating (I9.12j) from the hyperplane t = along the 
lines with t + r and cj = x/\x\ fixed: 

/•t+r 

(9.18) IzV^rw)! < / |dZ / /i 1 (t + r- / 9,pw)|aV>+ l^/i 1 (0, (t + r)w) | 



< Ce — — n ^ , ,, r^- + — — rr- ■ < 



r (1 + 1 + r) 1 -^! + |t - p|)i+7 (1 + t + r) 1 -^-* - (1 + 1 + ry- s (l + \t- r|)T 
A similar argument yields (|9.13p for r < t. Inequality (|9.14p follows from (|9.13|) using that \df\ < 

E,/|=i l^/l/(i + t + ' □ 

The next subsection is devoted to the preparational estimates for the inhomogeneous terms F^ v 
and F® u = D g h^ u arising in the equation (|9.5p for the tensor h l . 
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9.1 Estimates for the inhomogeneous terms 

These estimates will play a key role in the derivation of the improved decay and energy estimates in 
the following two sections. As we have mentioned above the quadratic terms in F^ u in (|9.ip do not 
satisfy the standard null condition. Nevertheless, a special tensorial structure and the wave coordinate 
condition will allow us to obtain favorable estimates. 

Recall that the inhomogeneous term F^ v has the following structure: 

(9.19) F^{h){dh,dh) = p{d^d v h)^Q^ v {dh,dh) + G^{h){dh,dh), 

(9.20) p(d^, d v e) -.= \d^d P e aP - ^trvr d v tr9, 

The quadratic term Q[iv is a linear combination of the null-forms and G fll/ (h)(dh,dh) is a quadratic 
form in dh with the coefficients - a smooth function of h vanishing at h = 0. The following result is 
an immediate consequence of (|9.20p . see also Lemma 14.21 



Lemma 9.6. The quadratic form P satisfies the following pointwise estimate: 

(9.21) \P(d7r,d9)\Tu < \dir\\d9\ + \dn\\d0\, 

(9.22) \P(dn, < \dir \ru\de\ru + \dn \cc\d6\ + \dir | \dd\ CC 

Using the additional estimates on the hn component, derived in Proposition 18.21 under the as- 
sumption that the wave coordinate condition holds, we obtain the following: 

Corollary 9.7. Let metric g satisfy the wave coordinate condition (|8.1|) relative to coordinates {x^}, 
the quadratic form P obeys the following estimate on a 2-tensor = g^ u — m^: 

(9.23) \P(dh,dh)\ T u<\dh\\dh\, 

(9.24) \p(dh, dh)\ < \dh\ 2 TU + \Bh\ \dh\ + \h\ \dh\ 2 

In addition, assuming that \Z J h\ < C for all multi-indices \ J\ < \I\ and vector fields Z G Z, 

|z 7 p(a/i,5/i)| < Yl {\dz J h\ TU \dz K h\ru + \dz J h\\dz K h\) + \dz J h\ LT \dz K h\ 

\.J\ + \K\<\I\ \J\+\K\<\I\-\ 

+ \dZ J h\\dZ K h\+ Y \Z' h h\\dZ- h h\\dZ- h h\ 

\J\ + \K\<\l\-2 [Ji|+|Ja[+|J 3 |<|/| 

Proof. The estimate (|9723|) follows directly from (pOT|h To prove (|9T24"|) we use (f9~22|) and that by the 
wave coordinate condition \dh\LL ^5 \dh\ + \h\ \dh\. 

We now note that Z I P(d tl h, d v h) is a sum of terms of the form P{d a Z J h, dpZ K h) for some a, j3 
and \J\ + \K\ < I: 

|Z 7 P(5/i,a/i)| < C \P(dZ J h,dZ K h)\ 

\J\ + \K\<\I\ 
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It follows from ([02]) and Proposition O that 



(9.25) Yl \P(dZ J h,dZ K h)\< \dZ J h\ T u\dZ K h\ru + \dZ J h\ LL \dZ K h\ 
|J|+|K|<|/| |J|+|K|<|/| 

< Y \dZ J h\\Z K h\ + \dZ J h\ T u\dZ K h\ TU 

\J\ + \K\<\I\ 

+ E ( Ei^'^ + E \9z j "h\ + y \z- h h\\dz j -h\)\dz K h\ 

\J\+\K\<\I\ |J'|<|J|-1 |J"I<|J|"2 |Ji[+|J2|<|J| 

which proves the result. □ 

We now state the complete estimates for the inhomogeneous term F^y. 

Proposition 9.8. Let F^ v = F^ u (h)(dh,dh) be as in (|9.19p and assume that the wave coordinate 
condition (18.11) holds for the metric = m^ v + relative to coordinates {x^ 1 }. Then 

(9.26) \F\ TU < \dh\ \dh\ + \h\ \dh\ 2 , 

(9.27) \F\<\dh\ 2 TU + \dh\\dh\ + \h\\dh\ 2 

In addition, assuming that \Z J h\ < C for all multi-indices \ J\ < \I\ and vector fields Z £ Z, 

(9.28) \Z I F\< Y {\dZ J h\ T u\dZ K h\ T u + \dZ J h\\dZ K h\) + Y \dZ J h\ LT \dZ K h\ 

\J\ + \K\<\I\ \J\ + \K\<\I\-l 

+ Y \dZ J h\\dZ K h\+ Y \Z J:i h\\dZ J ' 2 h\\dZ Jl h\ 

|J|+|if|<|/|-2 \Ji\ + lh\ + \Js\<\I\ 

Proof. First 

< C Y \Z Ik h\- ■ ■\Z h h\\dZ l2 h\\dZ h h\. 

|/i|+...+|/ fe |<|/|,fe>3 

Since ZQ{du,dv) = Q{du,dZv)+Q{dZu,dv)+a^Q ij {du,dv), and \Q^(dh, 8k)\ < \dh\ \dk\ + \dk\ \dh\ 
it follows that 

|Z J Q M „(dM/t)| < C Y \Q^(dZ J h,dZ K h)\<C Y \dZ J h\\dZ K h\ 

|j|+|fe|<|/| |j|+|fc|<]/| 

and the proposition follows. □ 

The next step is to estimate the extra inhomogeneous term O g h° appearing in the wave equation 
O g hj lu = Fftv — Ogh°^ u via the decomposition h = h 1 + h°. 

Lemma 9.9. Let 

F° u = U g h%, h% = x (r)x(t/r)M 6^ 
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Then, if the weak decay estimates in Corollary '\9-4\ hold, we have 



| Z - F 0|< C^( 4+ I,, + + M ,-*, ,>0 _ 2 

|Ce(t + |g| + 1)" 3 , q<0, 

More generally, 

^o^j^ + M + l)- 4 , q >0, C N e E , 

Proof. Set 

F° = B J ,/ l o = F 00 + F 01 , ^^D/! , F m = H^d a d p h° 
One can easily check that 

IZ'F 00 ! < Cn£ .„ , and F 00 = 0, for r < t/2 or r > 3t/4 

1 1 - (t+ \q\ + l) 3 ' 7 7 

'-(t + M + i) 3 ,^, 1 

On the other hand, with the help of (|9.13p for h°, 

y \z j h\ < Cj f + y \z j ^h i \...\z j -h\ 

|J|<fc 11 |Ji|+...+|J n |<fc 



since H = —h + 0(h 2 ), h = h° + h . Using f)9.13f) . this time for h 1 , we obtain 

CW^Ct+lgl + ir^Cl + lgDVa, g < 0, 



□ 



10 Decay estimates for the Einstein equations 

In this section we put to use the wave coordinate condition (j9.9p for the metric g(t) in coordinates 
{x /x } At= o,....3, which, as explained in the Remark 19,21 is satisfied on the maximum interval of existence 
[0,Tb], and the decay estimates for solutions of the wave equation Ugcf)^ = W^ u , derived in Corollary 
\7.2\ to upgrade pointwise estimates of Corollary 19,41 

Proposition 10.1 (Estimates for h). Let h = h 1 + h be a solution of the reduced Einstein equations 
(|9.ip . Assume that h 1 verifies the energy estimate (|9.10p on the time interval [0,T]. Then for any 

te[o,T], 



10.1) \dh\ LT + \dzh\ LL < 

10.2) \h\ LT + \Zh\ LL < 



Cs(l+t + \q\)~ 2+s (l + \q\)- s , q>0 

Ce{l+t + \q\)- 2+& {l + \q\) l l 2 , q<0 

'Ce(l + t + \q\)-\ q>0. 
Ce(l + t+\q\)- 1 (l + \q\)V 2 + 5 , q<0. 
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Furthermore, 

(10.3) \dh\ TU <Ce{l + t+\q\)-\ 

(10.4) \dh\<Csr x \nt. 

Proposition 10.2 (Estimates for h 1 ). Under the same assumptions let 7' < 7 — 5 and fjf > 8 > be 
fixed. Then there exist constants M^, and e^, depending on ("f',u',5), such that 

(10.6) | Z V|<( C ^ 1 + ( + l9l) tr (1 + kl) :I:. 5> ° , \I\=k<N,2 + 2 
The same estimates hold for h° if we replace 7' by M^e. 

Remark 10.3. Note the difference between the estimates (|10.3|) . (|10.4p for h and f 1 1 . 5 1) with |J| = for 
h 1 . 

The estimates (11 . 1 [ >- f fl . 2 f) follow from the wave coordinate condition combined with the weak 
decay estimates of Corollary 19.41 In the derivation of the sharp decay estimates f)10.3j) - ()10.6j) we will 
use decay estimates of Corollary 17.21 with various weights. 

10.1 Proof of flHEUHjHEZ) . 

We rely on Proposition 18.21 to establish an even more general version of the desired estimates. We 
remind the reader that under the assumptions of Proposition 110. ll both the wave coordinate condition 
(|9,9[) for tensor g^ u and the weak decay estimates of Corollary 19.41 hold true. 

Lemma 10.4. Under the assumptions o f Proposition [10. 1\ 

(10.7) 

' e{l+t + \q\)' 2+25 (l + \q\)- 25 , q>0 
+ t + \q\)- 2+25 (l + \q\) 1/2 - 5 , q<0 



Y,\dz I h\ LL + £ \dZ J h\ CT < £ + 

\I\<k \J\<k-l \K\<k-2 



10.8) 

'e(l + t+\q\)-\ q>0 



Y^\z i h\ LL+ £ \z j h\ CT < E / \ 9zKfi \ 

m<* \j\<k-i \K\<k-2 s ^ =const 



e(i + t + | (? |)- 1 (i + kl) 1/2+5 , g<o 



Here the sums over k — 2 are absent if k < 1, the sums over k — 1 are absent if k = and J 

s,uj=const 

stands for an integral along the segment r + |y|,y/|y| = const connecting a given point (t,x) with a 
hyperplane t = 0. 

Proof. The proof follows immediately from Proposition 18.21 and the estimates of Corollary 19.41 since 
(10.9) 

Y / \dZ I H\ LL + E \ dzJ H\ LT < ^\dZ J H\+ E \9Z K H\+ £ \Z h H\\8Z h H\. 

\I\<k \J\<k-l \J\<\I\ \K\<k-2 h+h=I 

Integrating fjlO. Tj) along the lines on which the angle u = y/\y\ and the null coordinate s = t + \y\ are 
constant, as in the proof of Corollary 19.41 and using (|9.13j) at t = 0, yields (|10.8p . □ 
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10.2 Proof of (gEED - pa) . 

We apply the L°° estimates, derived in Corollary \7.2\ for the reduced wave equation 



with F^ v given in (|9.2|) or f|9. 19j) - f)9.20|) . Observe that the weak decay estimates of Corollary [9 
guarantee that the tensor H^ u = g^ v — m^ v verifies the required assumptions of Corollary 17.21 



First we derive the L°° estimates for F^ v . 



Lemma 10.5. Suppose that the assumptions of Proposition UU. 1\ hold and let F^ v = F fiu (h)(dh,dh) 
be as in (19.21). Then 



(10.10) \F\ru < Cet-V 2+S \dh\, 

(10.11) \F\ < Cet- 3 / 2+5 \dh\ + C\dh\ 2 TU 

Proof. This follows from Proposition 19.81 and the weak decay estimates of Corollary 19.41 □ 

The 2-tensor h is a solution of the wave equation Ogh^ = F^ u . We now recall that that according 
to Corollary 17.21 with a trivial weight w{q) = 1, i.e., a = 0, we have that for any U, V £ {L,L,A,B}, 



(10.12) (l + t+\x\)\\dh(t,x)\uv<Bvp 

°^|/|<i 

+ ! ^l + T)\\\F(T,-)\ uv \\ Lx{DT) +Y,^ + r)- 1 \\Z I h(T,-)\\ Lx{DT) yr. 

|/|<2 

Then with the help of Lemma 110.81 and the weak decay estimates of Corollary 19.41 we obtain 
Lemma 10.6. With a constant depending on 7 > we have 



(l + t)\\\dh\ TU (t,-)\\ L o <Ce + Cs / (l + T) s - 1 / 2 \\dh(T,-)\\ L o a dr, 

Jo 

(1 + t)\\dh(t, OHioc <Ce + cJ (e(l + T ) 6 -V 2 \\dh(T, + (1 + r)\\\dh\ TU {r, -)\\ 2 L J) dr. 

Here, by assumption 5 < 1/4 so 5 — 1/2 < —1/4. The estimates (110. 3ft and (jl0.4[) are now 
consequences of the above lemma and the following technical result applied to the functions b(t) := 
(1 + t)\\ \dh\ru(t, and c(t) : = (1 + t)\\dh(t, 

Lemma 10.7. Assume that the functions b(t) > and c(t) > satisfy 

(10.13) b(t) < Ce (f o (l + s)- 1 - a c(s)ds + l^j 

(10.14) c(t) ^ Ce (J o (l + sr 1 ' a c(s)ds + lj +C J (l + s)- 1 b 2 (s)ds 

for some positive constants such that a > C 2 e and a > 4Ce/(l — 2Ce). Then 

(10.15) b(t) < 2Ce, and c(t) < 2Ce(l + oln (1 + 1)) 
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Proof. Let tq be the largest time such that holds. Substituting these bounds into (|10.13|) - (|10.14p and 
taking into account that 



T ) e ~ T dr = 2a' 1 + 1 



/*oo /*oo 

/ (l + s)~ 1 - a (l + aln(l + s))ds<a^ 1 / (1 
Jo Jo 

we obtain that for any < t < tq 

b(t) < Ce(2Ce(l + 2a" 1 ) + 1) < 2Ce, 

c(t) < Ce(2Ce(l + 2a- 1 ) + C 2 eln(l + t)) < 2Ce(l + aln(l + t)), 
which implies that tq = oo, as desired. □ 

10.3 Proof of (flTO -f lTM!) 

The proof of the estimates for the tensor h 1 proceeds by induction. We assume that (I10,5|) -( fl0.6p hold 
for all values of multi-index \I\ < k and prove the estimate for |/| = k + 1. The argument below will 
also apply unconditionally to the base of the induction k = 0. 
Once again the first step is to establish L°° estimates for Z 1 F. 

Lemma 10.8. Suppose that the assumptions of Proposition UU. 1\ hold and let Fu U = F fiu (h)(dh,dh) 
be as in (19.21). Then 



(10.16) \Z I F\<Ce l -^f^+C Yl \dZ J h\\dZ K h\ 

\K\<\I\ + |J|+|Jf|<|J|,|J|<|Jf[<|/| 

Proof. The result follows from Proposition ES] with the help of (j9llZj) - (f97L4"l) and fTlQ73|) . □ 
Recall now that the 2-tensor h = g — m — is a solution of the reduced wave equation 

(10.17) O g hl v = F^ = F tlv -F^ u , where F^ v = U g h% 

and that the terms F^ v and F® v have been treated in Proposition 19.81 and Lemma 19.91 respectively. 

To prove the estimates for Z 1 ^ with vector fields Z G Z we commute the equation (jlO. lTj) with 
Z 1 . By Proposition 1531 



(10.18) p^/i 1 ! < li'i? 1 ! + (i + ty 1 Y E \Z J H\\dZ K h 1 \ 

\K\<\I\, |J| + (|X|-1) + <|/| 

(10.19) 

E ( E E i^Ur+ E \z j "h\)\ 9 z^\ 

11 |A'|<|/| |J|+(|A-|-1)+<|J| |J'|+(|i<r|-l) + <|/|-l |J"|+(|K|-l)+<|/|-2 

where (\K\ - 1) + = \K\ - 1, if \K\ > 1, and 0, if \K\ = 0. Using Lemma [TUH and (|10.5j) . which is 
inductively assumed to be true for |/| < k, we get 

(10.20) (l + Mr 1 E (\z j h\ ll + \z j 'h\ lt +\z j "h\) 

\J\<k, I J'\<k-1, 1 J"\<k-2 

< {C k e{l + t+\q\)~ 1+M ^{l + \q\)- l - M ^ + e {l + t+\q\)-\l+\q\)- 1 , q > 
~lC7 fc e(l + t+|?|)- 1+M * e (l + M)- 1 /^' + £ (i + i+ | g |)-i(i + | g |)-i/2+5 ) g<0 
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while contribution of the terms in (|1U.19|) coupled to the highest order term - dZ K h}, \K\ = |/|, 
according to (|10.2p amounts to 

(10.21) (1 + kl)- £ (\Z J B\ll + \H\ LT ) < ( £(1 + 1 + kr ! (1 + Hr !' q > ° 

Therefore splitting the sum in ()10. 19|) according to whether \K\ < \I\ or \K\ = \I\ and using the 

inductive assumption fjl0.5f> for \K\ < \I\ = k + 1, we obtain 

(10.22) 

iSi 1+t 1 1 wi+t+ir^i+ir 1 ^, ?<o 

Now by Lemma 19.91 



.^OKl^ + M + D-^d + kir 5 , q >o, m<N _ 2 

- V £ (t + |g| + l)- 3 , g<0, 



Set 



(10.23) n k+l (t) = (i + t) £ MgjazWftOllx-, ^) = (!!t!1!i/2-M' 9> ° n • 

where u' > 5 and 7' < 7 — (5. Then we have established that for |J| = k + 1: 

(10.24) ^p^ 1 ! < (1 + t)- 2 {en k+l (t) + e 2 C|(l + t) 2M ^ + Ce{l + t)- 1 / 2 ^') 
The weak decay estimates of Corollary 19.41 imply that 

07(g) Z (t, a?) < < W/1 , . lU/2h , . lU/2 _u' n < e{l + t) , a>0 

|e(l + t+|g|) 1+d (l + |g|) 1/2 (l + \q\) 1/2 * , g<0 

provided that a = min (// — 5, 7 — 5 — 7') > 0. The decay estimate proved in Corollary I7.2I therefore 
shows that for some constant C = C(k) we have the inequality 

(10.25) n k+1 {t) <Ce + C f\l + t)' 1 {en k+l (r) + e 2 (l + t) Ce + Ce(l + t)- 1 ' 2 -"') dr 

Jo 

The bound n k+ i(t) < 2Ce(l + t) follows from the Gronwall inequality. This proves (|10.5p . The 
estimate (j!0.6p then follows by integrating (I10.5P along the line uj = y/\y\,T + \y\ = const from the 
hyperplane t = using (|9,17p . 

11 Energy estimates for Einstein's equations 

Recall the definition of the weighted energy 

(11.1) £ N (t)= sup V f IdZWfwb), 

0<r<t mN M 
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where 

ni9 ^ Jl + (1 + M) 1+2 T, q>0 , r(l + 2 7 )(l+M) 2 ^ q>0 

(11.2) w = < „ , w = < , _ 

ll + (l + |g|)-^, g<0 ^^(l + lgl) 1 "^, g<0 



Recall also our decomposition 



gnv(t) = m^u + ^(i) = m^v + ^(*) + h^{t), h^t) = x( + )x(r) — 6^ 



of a local in time smooth solution guvif) of the reduced Einstein equations, and the definition of the 
tensor H» v = g^ v - . 

In this section we prove the following result. 

Theorem 11.1. Let g^ u (t) = h^ u (t) + m^ u be a local in time solution of the reduced Einstein equations 
(19, 1[) ) satisfying the wave coordinate condition (|9.9p on the interval [0,T). Suppose also that for some 
< // < 1/2 and < 7 < 1/2 we have the following estimates for < t < T, all multi-indices 
\I\ < N/2 + 2 and the collections T = {L, S\, S%}, U = {L,L, S\, S2}: 

(11.3) \8H\tu + (1 + Wir'lHhc + (1 + Wir'lZHlcc < Ce(l + t)~\ 

(n.4) + ^ + 1 + * t i g| \Bz<h\ < ( cg(i + * + |g|r :^? + |9i r?; 2 /, q>0 

l + kl l + kl ~\C£(l + i + |g|)- 1 + cre (l + |g|)- 1 /2+M' J g<0 

(11.5) £jv(0) + M 2 < e 2 . 

Then there is a positive constant c independent of T such that if e < c~ 2 we have the energy estimate 

(11.6) £ N (t)<C N e 2 (l + t) C£ , 
for < t < T. Here CV is a constant that depends only on N . 



Assuming the conclusions of Theorem 111.11 for a moment we finish the proof of the Main Theorem 



II. 1 End of the proof of Theorem 19.11 

Recall that T was defined as the maximal time with the property that the bound 

£ N (t) < 2C N e(l + t) s 

holds for all < t < T. Assuming the energy bound above we have establsihed in Propositions 110.11 
110.21 the decay estimates for the tensors h = h + h l and h 1 respectively. Direct check shows that the 
estimates of Propositions I10.T1I10.2I imply the assumptions (lll.3p - (111.4|) . The conclusion of Theorem 

III. 11 states that the energy 

£ N (t) <C N e 2 (l + t) C£ , V0<t<T 

Thus choosing a sufficiently small e > we can show that £jv(i) < CVe(l + t) s thus contracting the 
maximality of T and consequently proving that g^ is a global solution. It therefore remains to prove 
Theorem 111.11 
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11.2 Proof of Theorem [TTTT1 

Proof. Recall that the components of the tensor h^ u = g^ u — satisfy the wave equations: 

(11.7) D g h liU = F lw „ 

Fpv = P(d„h, d v h) + Q^(dh, dh) + G^ v {h)(dh, dh), 

(11.8) Pid^h, 8 u h) = ^m^m^'d^hap d u h a , p , - \m aa ' m^'d v h m . 

Our goal is to compute the energy norms of Z h 1 , where Z G Z and h 1 is a solution of the problem 

(11.9) U g hl v = Fl v , where F 1 = F — F°, F° = a g h° 
Commuting with the vector fields Z 1 , Z = Z + cz we get 

(11.10) UgZ'h^^Fll, 
where 

(11.11) F 11 = Z 1 F — Z 1 F° — D 1 , D 1 = [Z I U g h l -U^h 1 ) 

We base our argument on the energy estimate (16. 9j) for a solution of the wave equation \D g cp = F 
of Proposition 16.21 Observe that the conditions of our Theorem on the tensor H^ v = g^ v — m fiv imply 
the that the assumptions of Proposition 16.21 for the metric g hold true. In particular, we have 

(11.12) f \dcf>\ 2 w + f f \d<j>\ 2 w'<&! \d<j>\ 2 w + wf I C^fJM! + p g( j>\ \d<j>\] w 
7s t Jo Jt, t Jt, Jo Jz t v l + * ' 

This applied to (jll.lOp gives 
(11.13) 

! \dZ I h 1 \ 2 w + f / IdZWfw'KS [ Idh^w + W [* f ( C ^ Zlk ^ + \F(\\dZ I h 1 \)w 
Jz t Jo J s r Jy, Jo Jz t ^ 1 + * ' 

< 8 [ \dh 1 \ 2 w + 16 f [ ( — \ QZlhX \ 2 w + £ -i(|Z J F| 2 + \D J \ 2 ) (l + t)w+ |Z J F°| {dZ'h'l 



We begin with the following estimate on the inhomogeneous term F: 
Lemma 11.2. Under the assumptions of Theorem \11.1\ 

mu) \z i f\< V ( £ 1^M i gO + MT^! \- B7 J h | | ^ 2 l^il 
(1L14) |Z Fl ~ ^ ^~TT^ + (TTi+klF^ |9Z 7111 + ITiTg TTM 

^ e|9Z J /n| e 2 



|J|<|J|- 



^{1 + ty-Ce (l + t+\q\r 
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Proof. According to Proposition 

\Z'F\< Yl {\dZ J h\ T u\dZ K h\ TU + \dZ J h\\dZ K h\)+ Yl \9Z J h\\dZ K h\ 
\j\+\k\<\i\ |j|+|i<r|<|7|-i 

(11.15) + \Z j3 h\\dZ j2 h\\dZ Jl h\ 
|Ji|+|J 2 |+|J 3 |<|/| 

Since h = h 1 + h° and h° obeys the estimates \dZ J h°\ < e(l + t + \q\)~ 2 and \Z I h°\ < e(l + t + 
for all |/| < N, it follows that both h° and h 1 obey the estimates (|11.3|) and (jll.4p of the theorem. 
Substituting h = h° + h 1 into (|ll,15p and expanding the products we obtain terms in which either 
all factors contain h°, in which case we can estimate them by e 2 (l + 1 + |g|) -4 , or at least one factor 
contains h 1 , in which case we can simply estimate the other factors by (111 .3[ >- f Tl 1 .4f> . It follows that 
f 1 1 1 . 1 5 [) leads to the estimate 

j v , e\dZW\ ejl + W-^ - j x e 2 \Z^h^\ 

el^Z- 7 /! 1 ! e 2 



EC /(, i 



□ 



Lemma 11.3. Under the assumptions of Theorem \ll.l\ 

e- 1 f hz'Ffil + QwdxdtZ Y f j e ^f^ 2 w + \dZ J h l \ 2 w')dxdt 
Jo J \J\<\i\ + ' 



^ r T r \dz J h l \ 2 J J 3 

+ /„ / £ ( 1 + t y~2Ce wdxdt + e . 



\J\<\I\-l 
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Proof. The estimate is a straightforward application of Lemma 111.21 We took into account that 
w < w'{\ + \q\)(l + q^-) 2 ^ and the inequality \i < 1 — 2/x'. The estimate 

i \z J h x \ 2 , , /■ laz- 7 /! 1 ! 2 

■wax 5- / -7 T-TT to ux 



i + t + kl (i + kl) 2 ~ J (i + t + 

is the Hardy type inequality established in Corollary 113.31 of Appendix B. □ 

Next we estimate F° = \J g h°: 
Lemma 11.4. The following inequality holds true: 

\dZ J h>\ 2 J ^ f T f f laryJull2 J \V2 eft 



jf y |Z 7 F°| laZ^V^i < Cjv£ ^ (j[ / wdxdt+ J o ( J \dZ J h l \ 2 wdx^j 

where Cjv denotes a constant that depends only on N. 
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Proof. We start with the estimate 

£ J {Z'F^ldZ^lwdxdtK, J 7 [J |Z J F°| 2 u;dx) 5 (J \dZ* h l \ 2 w dx) * dt 



Now by Lemma 19.91 

|Z / F 0| < fCiv£ 2 (t+kl + l)- 4 , q>0, C N e j t 

" \c^ £ (t+M + l)- 3 , g<0, (* + ^ + 1 ) 3 |^f/| 

Therefore, once again using a Hardy type inequality of Corollary 113.31 we obtain 

1*^0.2 j f f e 2 r 2 dr „ f°° e 4 r 2 (l + r) 1+2 ^ dr C N e 2 f \Z l h x \ 2 , 



< C ^ g2 + 



r2 

(1 + t)" 3 ' (1+t) 4 

and the result follows. □ 



J \dZ I h 1 \ 2 wdx 



The estimate for the term containing the commutator term D = OgZ^h} — Z I \Ji g h 1 follows from 
the following 

Lemma 11.5. Under the assumptions of Theorem \11.1\ 
(11.16) e" 1 j J \D g Z I h 1 -Z I n g h 1 \ 2 (l + t)wdxdtii 

^ e E f f(^^™ + \dZ J h^ W >)dxdt + e ^ f fjS^vdxdt + e* 
\J\<\i\ |J|<|/|-1 ^ ' 

We postpone the proof of Lemma 111.51 for a moment and finish the proof of Theorem 111.11 
Using ([11.13P together with Lemmas lll.3Hll.5l yields 

(11.17) 

+ Ce E / T / ( 'df^ 1 ^ + ISZ^p^dxdi + Ce f I j^^^wdxdt + Ce 3 

\J\<\i\ + ' |J|<|/|-1 ' + ' 

where Cat depends only on |/| < N. As before we denote 

]2 w dx 



S k (t) = sup V / \dZ z h x \ 

0<r<t _ c , |r| JS T 



' zez,\i\<k' 

and let 



5 fe (t):= ^ f [ \BZ I h 1 \ 2 w'dx 



zez, \i\<k 



o js t 
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It therefore follows that 
(11.18) 

*M + m < SE t{0) + Ctf, W + £ °gfc>* + jf g ffifff <tr + C E ' + jf fftiW <ir 

For Ce sufficiently small we can absorb the space-time integral S k (t) into the one on the left hand- 
side at the expense of at most doubling all the constants on the right hand-side. Similarly, since 
16C N el/4Ejf 2 < E k + U 2 C 2 N e 2 and E k {t) is increasing, we can absorb 1/4 jj E k (r) (1 + r)~ 3 / 2 < 
\/2E k (t). If we also use the assumption -Ejv(O) < e 2 we obtain for e > sufficiently small 

(11.19) E k (t) + fc(t) < ^e 2 + f ^Mdr + f gg^ ^r 

where the last term is absent if A; = and Cat is a constant that depending only iV. 
For k = this yields the estimate 

^(t)<C^+/*^§Mdr 
Jo 1 + r 

and the Gronwall inequality gives the bound 

E (t) < c N (i+tr £ 

which prove .6[) for k = 0. 

Assuming (111 .61) for k replaced by k — 1 we get from (111.19ft 



(11.20) E k (t) < C N e 2 + I* ^Idr + f ? ^ 

K ' kK! ~ Jo 1 + r 7o (1 + r) 1 "^ 



which leads to the bound 

E k (t)<C N e 2 (l + t) 2c£ 

This concludes the induction and the proof of the theorem. □ 

11.3 Proof of Lemma 111.51 

Proof. Define the tensor 

:= H»" - Hg», = - X {- t )x{r)^ lv 

Observe that H$ coincides with the tensor —h°. Define also the wave operator □]. = □ + H^dadp. 
According to (|5.9p of Proposition 15.31 



(11.21) \n lZ 'k>-z'n lh i\< Y. E (ir^ + ^f) 1 ^'' 11 

\K\<\I\ |J|+(|iC|-l) + <|/| m m 



+ E( E + E ^> z,<h 

\K\<\I\ |J|+([Jf[-l) + <|J|-l m |J|+([K[-l)+<|/|-2 m 
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Our goal is to obtain the estimate for the quantity 

rT 



\J\<N 

Let us first deal with the terms in (jll.21|) with \K\ < N/2 + 1. In this case we use the decay estimates 

(E3D 

\dZ K h\ < Ce{\ + t + \q\)- l+C£ (l + \q\y 1/2+f1 ' < Ce(l + t + \q\y l+C£ (l + M)"^ 

guaranteed by the assumptions of Theorem 111 .11 provided that fi < 1/2 — //. It is clear that in this 
case it suffices to consider the expression 

(11.22) F [( \ ZJHl? v \Z J HAIc + \Z J 'hAIt+\ zKh i\ 2 \ e 2 (l + \q\)- 2C£ w dxdt 



(l + t+k|) 2 (1 + kl) 2 J {l + t+\q\Y-^e(l + q _)2^ 

^ . T r \Z J H X \ 2 e 2 wdxdt f T f \Z K H 1 \ 2 e 2 wdxdt 



rvj 



(1 + M) 2 1 + t J J (l+|g|)2(l + t )l-2C £ 

T r \Z J H 1 \ 2 CC + \Z J ' e 2 (l + \g\)- 2Ce w dxdt 



+ 



(1+M) 2 (l + t + \q\y- 2 ^(l + q^ 



with \J\ < k, \J'\ < k — 1, | if | < k — 2, where k = |/|, After applying the Hardy type inequalities of 
Corollary 113.31 the above expression is bounded by 



e 2 



J / ^ Z i + t wdxdt + J / (f + t y-lc s w dxdt + J j {\dZ J Hi\ 2 cc + \dZ J> H^l^wdxdt 



where 



w 



w = mm [w , 



(i + t + \ q \y~ 2 ^'- 

Ignoring the difference, which we shall comment on at the end of the proof, between the tensors 
= g^ u — — Hq V and h 1 ^ = g^ u — m u ^ — h®, we see that the first two terms are as claimed in 
the statement of the lemma. We now recall that according to Lemma 115.41 of Appendix D 

Ys\ dZ ' ,H i\cc+ E \dZ J H 1 \ cr <Y,\dZ J H 1 \+ £ \dZ J 'H x \ 

\J\<k \J\<k-l \J\<k \J'\<k-l 

\Z J H X \ 



+ £ 



V ( \dZ J Hi\ + J 

l + i+M^V 11 ^1 



\j\< k - + t + ^ 



+ V 17* ff \\87^TJ I i Cg Xo(V2 < r/t < 3/4) Ce 2 

It is clear that in the sum Sl^l+ijal^ above at least one of the indices is < A^/2 and therefore we can 
use the decay estimates (fTPj) : for \J'\ < N/2 + 2 

+ \dZ J 'h\ < Ce(l + t + k|)- 1+Ce (l + \q\y 1/2+ »' < Ce(l + t + Mr 1+Ce (l + \q\)~ Ce '■ 

1 + \q\ 
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Thus, with \J\ < k, \J'\ < k - 1, 



(11.23) jf J (\dZ J Hi\ 2 cc + \dZ J 'H 1 \ 2 CT jwdxdt< J ^ \dZ J Hi\ 2 w' dxdt 



-e 2 




\j\<k 



^ f T f {dZJ'H^ 2 f T f f Ce 2 X0 (l/2<r/t<3/4) Ce A x 

,JtJ° JWTW^ wdxdt + Jo Jy (T+t+W (iTRTW^ 

where xo(l/2 < r/t < 3/4) is the characteristic function of the set where t/2 < r < 3i/4. Using the 
properties of the function w' , in particular that w 1 < «V(1 + \q\) we obtain that the above has a bound 

^ f T f \dZ J 'H x \ 2 J ^ f T f Ce 2 X 2 (l/2<r/t<3/4) J . 
+ > / / -J N1 ' wdxdt+ / / uv / ; , t J ' dxdt 

^ L J (l + t )l-2C7 e y Q _/ g)4(l + g |)l+a^ 



IJ'l<fc-l 



T r 4 

^(l + |g|) 1+27 dx^ 



+ + 

Once again we use the Hardy type inequality of Corollary I13.3I with a = to conclude that 

\J\<k 



J J (idZ'HtHc + ldZ^tllrJw'dxdt^ J Q J (\9Z J Hi\ 2 v/ + \+l w)dxdt 

T f _\dZ J 'Hx\ 2 

\J'\<k-l ' " 



as desired, modulo the difference between fli and /i 1 , in Lemma [11.51 To finish the proof of the Lemma 
for the case \K\ < N/2 + 1 it remains address the difference between the tensors 

= - + h ofM/ and h\ v = 9iiV - - h V 

The above expressions imply that 

(m + h° + h%l = (m-h° + H x y v 

Therefore, B{ v = -h lfiu + 0^ u ((h° + h 1 ) 2 ) and it follows that 



\dz J R x \ < \dz J h x \ + v (jz^Hdz^/i 1 ! + (z^^pz- 72 /! 1 ! + laz- 71 /! !!^- 72 ^ 1 ! + \z Ji h°\\dz j2 h°\), 

iai+ij 2 |<iji 

\dz J H x \ < \dz J h x \ + (\z Ji h°\\dz j2 h l \ + {z^h^idz^h^ + 1 a^- 71 ^° 1 1 /i 1 1 + |z Ji /i°||dz j2 /i°|) 



|Jl|+|J2|<|J| 
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Taking into account that for % = 0, 1 and | J'\ < N/2 + 2 

J^-j^ + \dz J 'h\ < ce(i + i + M) -1+Ce (i + \q\r 1/2+f1 ' < ce(i + 1 + kir 1+Ce (i + M)~ Ce - 

we obtain that 

Thus, with the help of the inequality w' < w/(l + \q\), 
|«7|<fe,[^|<fe— a + ^ + 



\A< k \ K \< k ~^- 



T <■ 1 \Z J h l \ 2 J Jm f T f 1 \Z K h l \ 2 , u , _ 2 



U y (i+t)(i + ki)^^ + y y (i +t) i-2c e( i +k |)2«'^+^ 

< ^ / / - ! wdxdt+ ^ Z ^_l Ce ^jwdxdt+ f f \dZ J h 1 ^ w' dxdt + e 2 , 

where to pass to the last inequality we once again used the Hardy type inequality of Corollary 113.31 

Returning to (|11.2ip we now deal with the case \K\ > N/2, which implies that \J\ < N/2 + 1 and 
allows us to use the decay estimates (jll.3p -( jll.4j) for Hi = —h 1 + 0(h 2 ). Therefore, the contribution 
of the terms with \K\ > N/2 to (Qi^/i 1 - PUih l \ can be bounded by 

\Z J H X \ , ^H^cc + WcT^xrKuU , V- l^llm^li 



y y( i z gjL + If Hllcc + MllcT )\dz K h'\+ y l -^4\dz 



(l + t) 1 -^ 

jc|=|j| |^|<|/| v 7 



E ^+ E 

fC]=|J' 

and the desired result follows. 

To estimate the commutator D^Z^/i 1 — Z I \Ji g h 1 it remains to address the term 

iH^dadaZ'h 1 - Z l (Bfd a daZ I h 1 \\<e — ^ V c) 82 ^- 

1 P V " P Jl ~ l+t + \q\ ^ 1 + kl 

Therefore 

[ T [ IHfdadaZ 1 ^ - Z^HfdadgZ^^il + ^wdxdt <e 2 y f [ Igg^f \ w dxdt 
Jo J \j\<\j\ Jo J 1 + t 



□ 
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12 Appendix A. Commutators 

Recall the family of vector fields Z = {d a , Q, a p = —x a db + xpd a , S = tdt + rd r }. In this section we 
address commutation properties of the family Z with various differential structures. Recall that for 
any Z 6 Z we have [Z, □] = — c^D, where cz is different from zero only for the scaling vector field S, 
c s = 2. 

Lemma 12.1. Let Z £ Z and let the constants c« be defined by 

[d a ,Z] = cJ l d fl , cJ? = d a Z» 
Then cll = c— = 0. In addition, if Q is a null form, then 

(12.1) ZQ(d<f>, dip) = Q(d</>, dZifj) + Q(dZcf>, dip) + Q(d<f>, dip) 
for some null form Q on the right hand-side. 

Proof. Since Z = Z a d a is a Killing or conformally Killing vector field we have 

(12.2) d a Z p + dpZ a = fm a p 

where Z a = m a nZ" . In fact, for the vector fields above, / = unless Z = S in which case / = 2. In 
particular, 

L a L /3 d a Z fS = 0. 

If Ca is as defined above and c a p = c^m^p = d a Zp the above simply means that cll = c— = 0. 
which proves the first part of the lemma. To verify (112. ip we first consider the null form Q = Q a p We 
have 

zQ^idcj), dip) = Q a p(dz^, dip) + Q a p{d<p, dzip) 

+ [z, d a ]<pd P i> - a <f>[z, d a ]ip + [z, dp\4>d a ip - d a <p[z, d p \ip 

= Q aj3 {dZ<p, d^) + Q a p{d(p, dZip) - c^d^dpip - dpcpd^ip) - c^{d^d a ip - d^d^) 
= Q a p(dZ0, dip) + Q Q/3 (d0, dZi,) - c£Qtf(d<f>, dip) - c^Q^ipcP, dip) 

The calculation for the null form Qo(8<p, dip) = m a ^ d a (pdpip is similar and we leave it to the reader. □ 

For any symmetric 2-tensor ir and a vector field Z S Z define 

(12.3) 7rfd a dp = Tr a P[d a dp,Z], i.e., n tf - «<*r, f> -l^ <* 



"7 1 7 



Lemma 12.2. The tensor coefficients verify the following estimate: 

(12.4) \tt z \ ll < 2\ir\ LT . 
In general, 

(12.5) [^d a d p ,Z I ]= Yl K^dadpZ 1 *, 

h+h=I,\h\<\I\ 

\K\<\J\ K+Z=J \K\<\J\-2 

for some constants c^f^ and d"^^ u . Here the sum (|12.5|) means the sum over all possible order pre- 
serving partitions of the multi-index I into multi-indices I\, I^. 
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Proof. First observe that since the vector fields Z are linear in t and x we have 

2 

*/3j 



[<9^3, Z] = [dp, Z]d a + [d a , Z]dp = c^d^da + c^&ydp, 



which proves the first statement, while the second follows since cjj = 0. 
To prove fll2.5f) we first write 

Z I {^d a 8^) = {Z K ^)Z J {d a d^) 

K+J=I 

Then we observe that 

(12.6) Z J d a d^ = ]T [Z«i , [Z l h ,[..., [Z^-i , [Z^ , d 2 a0 \ ]...]] ] Z J ^, 

Ji+J'2=J, Ji=(ie 1 ,...,ii n ) 

where the sum is over all order preserving partitions of (1, k) into two ordered sequences (£\, £ n ) 
and (£ n +i, ■■■Ak) such that J<i = (te n+1 , t£ k ). It therefore follows that 

n Ja ? = - {Z K ^) I 2 *' [ Z " 2 > \- ZH -^ \ ZH i Q U\ -] ] ] 

K+L=J,L=(i, 1 ,...,i, l ) 

The desired representation follows after taking into account that 

(Z K ^)[Z,d 2 ap ] = -(Z K n a /)d a dp 

□ 

For a symmetric 2-tensor H and a vector field Z £ Z we set Z = Z + cz , where cz is the constant 
in the commutator [Z, □] = — c^D, and 

(12.7) H J ^= c[f^Z M H^ = -Z J H a ? - Z M Hf+ ^ d^Z M H" v 

\M\<\J\ M+Z=J |M|<|J|-2 

Corollary 12.3. Let U g = □ + H af3 d a dp. Then 



(12.8) U g Z4> - ZO g <P = -(ZH a ? + Hf)d a „pv, 

'A \ZH\ L _ 
+ t + \q\ l + \q\ 



(i2.9) p,z, - zu,,\ < ( Hf + \**^ + Wit , E ^ 



In general, 





\n g z<f>- zD g 


» 


- z'ngcfr = - 


\UgZU 


>-z I o g< f>\ < 


E ■ 


( E 



I/I<1 



Y H^dadpZ 1 ^, 

h+h=iAh\<\i\ 

it^m E E i^i 

|-K1<|i1, |J|+(|ic|-i)+<|/| 

+TTTa\ E ( E 1^1-+ E I^Ur+ E l^l)l^ 

|Jf|<|/| |J|+(|*:|-1) + <|/| l-^'l+df |-1)+<|/|-1 |J"|+(|K|-l)+<|/|-2 

w/iere (\K\ - 1)+ = |K| - 1 if \K\ > 1 and {\K\ - 1)+ = if \K\ = 0. 
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Proof. First observe that 

ZU 9 4> = {Z + Cz )D<f> +{Z + c z )H a ^d 2 a ^ 

= UZ<j> + H<*Pdl p Zcf> + (ZH^)d 2 a ^ + (Hf + c z H^)dl^ 
= D g Zcj> + (ZH^)d 2 a ^ + (Hf + c z H^)d 2 ap( p 

Recall now that the constant cz is different from only in the case of the scaling vector field S. 
Moreover, in that case 

Hf + c s li afi = 

The inequality ([12.91) now follows from (112. 8h . (j!2.4f) and the estimate (15. 1\ . The general commutation 
formula (|12.10|) follows from the following calculation, similar to the one in Lemma 1 12.21 We have 

Z 1 ^ = Z 1 ^ + ^H^d 2 ^ = uz 1 ^ + z J H^z K dl p( j ) 

J+K=I 

If we now use (|12.6p we get (jl2. 1Q|) as in the proof of Lemma 112.21 The inequality ([12. lip now follows 
from (112.101) . (THED and the estimate ([121). n 



13 Appendix B. Hardy type inequality 

In this section we prove a version of the classical three dimensional Hardy inequality 

2 



f ^-^-dxKA [ \Vf(x)\ 2 dx. 
7r3 N Jr 3 



The Hardy inequality converts the weighted L? norm of the function / into the L 2 norm of its gradient. 
These type of estimates prove to be useful in the context of energy estimates for solutions of a quasilinear 
wave equation Ug^cp = F, where the energy norms contain only the derivatives of (f> while the error 
terms, generated by the metric g(4>) also depend on the solution 4> itself. The disadvantage of the 
classical Hardy inequality in this context is that it requires a costly weight r~ 2 . The "cost" here refers 
to the rate of decay of the weight in the wave zone r ~ t. Therefore we seek a version of the Hardy 
inequality with the weight dependent on the distance to the cone r = t rather than the origin r = 0. 

Lemma 13.1. Let 0<a<2, 1 + [i > and 7 > 0. Then for any function u G Cg([0,oo)) and an 
arbitrary t > there is a constant C, depending on a lower bound for 7 > and 1 + fi > 0, such that 

f* u 2 r 2 dr f°° u 2 r 2 dr 

(13.1) / : + 



(1 + \r - t\) 2+ ^ (1 + t + r) a J t (1 + \r - t\y-"i (1 + t + r 



<C - ' , — , x —dr + C \d r u\ 2 ^—^ - s r 2 dr 



(1 + \r-t\Y {l + t + r) a J t 1,1 (l+t + r) a 

Remark 13.2. The inequality ()13.ip with a = appears to be the precise analogue of the classical 
(spherically-symmetric) Hardy inequality. The presence of the additional growing weight (l + |i — r|) 1+7 
seems to be necessary and in fact fits perfectly in the context of our weighted energy estimates. In the 
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previous works |L-R2j the estimates used to convert a weighted norm of a solution into a norm of its 
derivative, e.g., 



f (1+ i;_ t i)^ ^sc((i +t )V(i + Qi 2 + / 

were more reminiscent of the classical Poincare inequality. 
Proof. Set 

m(q) 

Then 

m'(g) = 

Since a < 2 



(l + |r-t|) 2 +^ 



+ g>0 
" \(1 - g)- 1 ^, q<0 

7 (l + g)- 1+ T, g>0 
(l + //)(l-g)- 2 -^, g<0 



V v ' v V Vr 1 + t + r mr-t)J 



a m'(r — t)\ r m{r — t) r m'(r-t) 



> 



-r 1 + t + r m(r - t) J (1 + t + r) a ~ (1 + t + r) a 
Hence 

d r (V 2 (l + i + r)~ a m(r - i)0 2 ) > m'(r - t)r 2 (l + t + r)~ a (p 2 + 2r 2 (l + f + r)~ a m{r - t)(pd r 4> 
If we integrate the above from to oo and use that <j) has compact support we see that 

' u \ 2 r 2 dr f 00 . , JO , r 2 dr 

m [r — t)4> t r— <2 I m(r — t)(po r <p- 



{l + t + r) a ~ Jo y ' ^(l + t + r) a 

Since m' > and m > it follows from Cauchy-Schwarz inequality that 

f°° .. xj2 r 2 dr r- f°° m(r-i) 2 ,., r 2 dr 

Jo (l+t + r) a Jo m'(r-t) (l + t + r) a 

from which the lemma follows. □ 
Corollary 13.3. Let 7 > and /i > and set, for q = r — t, 

'l + (l + kl) 1+27 , q>0, 



w(q) 



i + (i + kl)" 2/i , q<o 



Then for any — 1 < a < 1 and any <j> € Cj 



|^| 2 n;dx f \nj.\2 wdx 



{l + \q\) 2 (l + t+\q\y-°~J 1 ri (l + t+| g |)i-« 
i/ m addition a < 2 min (7, /x), 



H 2 (i + kl)~ a ^ /- mjLl2 . , w w 

' 11 1 mm ( it; , — : , , n _ a J dx 



(l + \q\) 2 (l + t+\q\y-* (1+q^ - J 1 Yl ' ' (l + t + 

where q- = \q\, when q < and g_ = 0, when q > 0. 
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14 Appendix C. Weighted Klainerman-Sobolev inequalities 

In this section we provide a straightforward generalization of the Klainerman-Sobolev inequalities, 
expressing pointwise decay in terms of the bounds on L 2 norms involving vector field Z G Z. We 
consider energy norms with the following weight function 



w = w(q) = 
for some < 7 < 1. Therefore, 

w ' =: w'(q) = 



1 + (1 + |g|) 1+27 , when q > 
1 + (I + \q\)- 2 ^ when q<0 



J (l + 2 7 )(l + |g|) 2 T, when ^>0 
[2fj.(l + M) _1_2m , when q<0 
w' < 4w(l + Iqiy 1 < 16 7 ~V(1 + q^. 



We have the following global Sobolev inequality 

Proposition 14.1. For any function (f> £ Cq°(IR 3 ) and an arbitrary (t,x), 

\<l>{t, x)\(l+t+ \q\)[(l + \q\)w(q)] 1/2 < C ^ Ww 1 ' 2 Z 1 <j>{t, .)\\ L *, q = t-r 

M1<3 

Proof. Fist note that it is sufficient to consider two cases when the support of <p is in the sets r <t/2 
and t/4 < r respectively. We argue as follows. Let x( r ) be a smooth cut-off function such x( T ) = 1 
when r < —3/5, andx{r) = for r > —1/3. Define 

Mt,x) = X { r -^Mt,x), fc(t,x) = (1 - xC-^))4>(t,x) 
r + t r + t 

The supports of functions 4>\,4>2 then belong to the desired regions and S|/|<3 \Z I( t>i\ < CS|/|<3 l-^"Vl> 
since \Z I \i\ < C- 

First for r < i/2 consider the rescaled function c<j>t(x) = <f>(t,tx), where \x\ < 1/2 in the support. 
By the standard Klainerman-Sobolev estimate 

ii&(z)Hl- < c <cY1 w(t^d a mM\\ L i 

\a\<2 \a\<2 



Since 



we obtain that 



|t-r|H|^0(t,x)|<C^ 

\a\<k \I\<k 



\m,x)\\ LT <ct-^ 2 Y J \\z I ^x)\\ Li . 
\i\<k 

Multiplying both sides of the above inequality by [w (— t)] 1 ^ 2 , which is ~ [ - w(?)] 1 ^ 2 when r < t/2, proves 
the proposition in the case when <f> is supported in the region r < t/2. 
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When r > i/4 we can estimate the L°° norm by the L 1 norms of derivatives: 
\<t>(t,x)\\l+t+\q\) 2 w(q)(l + \q\) <C Y [ [ \dZd q (w(q)(l + \q\)(t + q) 2 <f>[t, (i + q)uf) I dqdS(u>) 

Since + |g|) < Cw(q) and |g| < C\t + q\ = Cr on the support of <f> it follows that we have the 

bound 



|a|<2,fc=0,l 
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w{q){t + qf {qd q ) k d^4>{t,{t + q)u) dq dS '(w) < C w(g)|Z J <H*, x)| 2 cte 



|J|<3 



R3 



□ 



15 Appendix D. Wave coordinate condition 

In this section we provide the details on the estimates following from the wave coordinate condition 
(15.1) 9 M (^v4det^) =0 

for a Lorentzian metric g in a coordinate system {x M }^=o,...,3, leading up to the proof of Proposition 
18.21 We recall the definition of the tensor H^ v = g^ v — mf*" and below state the consequences of (|15.1|) 
in terms of estimates for H. 

We first observe that for any vector field X and a collection of our special vector fields Z € Z we 
have that 



;i5.2) 



z'd a x a = d a (z I x a + Y cj«z J x~<) =a a ( Y cj 7 Q ^ 7 ), 

\j\<\i\ ' \j<\i\ 



where cf^ are constants such that 



c j« = 5«, for \J\ = \I\ and 



I L 
JL 



0, for |J| = |/|-1 



The last identity is a consequence of the relation between cj J and the commutator constants c a p = 
[d a , Z]p for which we have established in Lemma 112.11 that cll = 0. It therefore follows from f)15. II) 
and (fT5~2|) that 



(15.3) 



H^:=Z I H lu/ +^2 c I J ^Z J H lv , with H^:=H t 
\A<\I\ 



fJbV 



tiH, 



satisfies 

(15.4) d^HW + Z I d lx O lxv (H 2 ) = 0, where 0^(H 2 ) = 0(\H\ 2 ). 
Lemma 15.1. Let H be a 2-tensor and let i?W be defined by (115. 3D . Then, for j = 0, 1; 

(15.5) ^|^Z 7 F| LL + Y \d 3 Z J H\ LT < Y, \&Z K H\+Y\d JHlI] \LT 

\I\<k |J|<fc-l \K\<k-2 \I\<k 
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Proof. It is easy to see that (|15.3|) implies that for any a 

^|crZ 7 F| LL < \d a Z J H\ LT + Y, \d a Z K H\+Y J \d a H [I] \ LL , 

\I\<k \J\<k-l \K\<k-2 \I\<k 

Y \d a Z J H\ Lr < Y \d a Z K H\+ Y \d a H [J] \ LT 

|J|<fc-l \K\<k-2 |J|<fe-l 

and (fT531) follows. □ 

Lemma 15.2. Let Z represent Minkowski Killing or conformally Killing vector fields from our family 
Z and assume the tensor H satisfies the wave coordinate condition (115.11) . Then for any multi-index I 
tensor , defined in f|15.3[) satisfies the following estimate 

\8H^\ LT < Y \dZ J H\+ Y \Z Ik H\---\Z h H\\dZ h H\. 

\J\<\I\ h + ...+I k =I,k>2 

Proof. The wave coordinate condition (18. ip can be written in the form 

8^ = 0, where & v = (rnT + H» u ) y/\detg\. 

It follows from (fT5T2|) that 

M E c I J ^Z J G^=0. 
\J\<W 

Decomposing relative to the null frame (L, L, A, B) we obtain 

d q { Y c I J ^Z J G lu )=d s [ Y, c]^Z J G, v )-A,d A { Y Cj^Z'G^). 
\J\<\i\ \J\<\i\ \J\<W 

We now contract the above identity with one of the tangential vector fields T v , T G {L, A, B} to obtain 

L-TT'd^G^+Y Cj^T 'd q Z J G, /v \< Y WG\ 
\A<V\ \A<V\ 

We now examine the expression 
L? T v Z J d q G lv = I? T u 8 q Z J (im lv + H JV )y/\detg\\ = Y Ll T ^ {{Z Jl H lu )Z- h yJ\detg\\ 

Jl + J2 = J 



since = U 1 T v m^ v = 0. The desired estimate now follows from the identity yj\ det g\ = 1 + f(H), 
which holds with a smooth function f(H) such that /(0) = and f(H) = -triT/2 + 0(H 2 ) □ 

We now summarize the above results in the following 

Lemma 15.3. Let g be a Lorentzian metric satisfying the wave coordinate condition (115. 1|) relative to 
a coordinate system {x^}^=o,...,3- Then the following estimates for the tensor H^ v = g^ v — and a 
multi-index I hold true under the assumption that \Z J H\ < C, for all \ J\ < \L\/2: 

(15.6) \()Z'H\r. T <( Y \VZ J H\+ Y \9Z J H\ + Y \Z h H\\dZ h H\ 

\J\<\I\ \J\<\I\~1 |/i|+||/ a |<|/| 

(15.7) \dZ J H\ LL < ( Y \9Z J H\ + Y \ dzJR \+ E \Z l2 H\\dZ h H\ 

\J\<\I\ \A<\I\~* |/i|+|Ia|<|J| 



56 



We now establish that the tensor 

(15.8) fff = - £C, H$" := - x (I) x ( r )^^ 

obtained by subtracting the " Schwarzschild part" Ho from H obeys similar structure estimates to those 
of H. 

Lemma 15.4. Let g be a Lorentzian metric satisfying the wave coordinate condition (|15.ip relative 
to a coordinate system {x M }^=o,...,3- For a given integer k > we have the following estimates for the 
tensor H^ u = g^ v — m^ u — Hq U , under the assumption that \Z J H\ < C, V| J| < k/2, 

(15.9) \ dzlR Acc+ S \ dzlRl \cr ~ \ BzIh i\ + \ dzlR i\ 

\I\<k |/|<fc-l |/|<fe \I\<k-2 



l+t+ <7 V l + t + 



191 |/|<* 



+ V \7 I H \ \f)7 J H I I Cg Xo(V2 < r/t < 3/4) 
+ ^ + ( i + t+ U|)2 + 7T 



Ce 2 



(l + t + |g|) 2 (l + |t| + |g|) 3 ' 

where xo(l/2 < r/t < 3/4) is i/ie characteristic function of the set where t/2 < r < 3i/4. 
Proof. We define the tensors 



H]n> — Zl Hfiv +Z C J^7^ J ^7 I/ ' ^A*" ~~ H/j,p ^-tvH, 

\J\<\i\ 

Ho\jtl = ^Hq^ +^ Cj ^ Z J H ^ u , H 0fiu = H 0flv y-tvHo 



\J\<\i\ 

Hi^l = Z 1 Hi^ u c jfJ Z J Hi lv , Hi^u = Hxpu y~%xHx. 

\J\<\i\ 

with the constants Cj 7 such that = 0, if |J| = |/| — 1. Observe that the tensor iJoJS is defined 
in such a way that 

d*H W=z J (d*H 0ia ,) 
Since the tensor satisfies the wave coordinate condition it follows that 

(15.10) d»HW+Z I d fl O^(H 2 ) = 0, where (H 2 ) = 0(\H\ 2 ). 

On the other hand, calculating using the definition of Hq, 

d»H 0fl „ = 2x'(r/t)x(r)M/t 2 5 l/0 

Therefore, since = H 0fll/ + Hi^ v and HjQ = H ^l + H^l we obtain 

d^H^l = -Z'd.O^iH 2 ) + 2Z I ( X '(r/t) X (r)M/t 2 )5 uQ 
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Arguing as in the proof of Lemma 115.21 we derive 



(15.11) J2\ dzlR i\cc+ E \ dzlH AcT^Y.\ dH nLT+ E \ dzlR ^\ 

\I\<k |J|<fc-l |J|<Jfe |/|<fe-2 

(15.12) < ^ + E E \Z ln H\ ■ ■ ■ \Z h H\ \dZ h H\ 

\I\<k \I\<kh+...+In=I,n>2 

f15 13l + V \87 I H I l CeXo(l/2 < r/t < 3/4) 

(15.13) + ^ \BZH X \ + (1 + t+k |)2 • 

[J|<fc-2 V WJ 

where xo(l/2 < r/t < 3/4) is the characteristic function of the set where t/2 < r < 3i/4. Using the 
condition \Z I H\ < 1, V|/| < k/2 and the estimates 

\Z I H \ + (l + t+ {qDldZ'Hol < 1+f g +|g| 

we obtain 

(15.14) Y \Z In H\ ■ ■ ■ \Z h H\ \dZ h H\ < ^ \Z J R\ \dZ K H\ 

h+...+I n =I, n>2 |J|+]i<r]<|J| 

< J2 \Z J H \\dZ K H 1 \ + \Z J H 1 \\dZ K H \ + \Z J H 1 \\dZ K H 1 \ + \Z J H \\dZ K H \ 

\J\+\k\<W 

\z J n x \ \ e 2 



, , • Ill V" -' ' l + t+\q\) ' (l + t + \q\f 

and the lemma follows. □ 

16 Appendix E: L 1 — L°° estimates and additional decay 

The bounds 

(16.1) £ N (t)= \\w ll2 dZ I h\t r )\\ L 2 + \\w 1 l 2 dZ I ^{t,-)\\ L ,<C N (l + t) c ^, 

\I\<N 

with 

'(1 + |r-t|) 1 /2+7 ) r>t 
1, r<t 
imply, with the help of the weighted Klainerman-Sobolev inequality, that 

(16.3) i^^i + i^^i < + * + r)_ ;:r (1 + " " r|r !r' r>t > i^*-* 

K ' 1 11 ^- \c' N e(l + t + r)- 1+c » e (l + \t-r\)- 1 / 2 , r <t 

The exponent in the interior decay estimate ()16,3p in \t — r\ can be improved from —1/2 to —1 using 
the estimate 

(16.4) IZ 7 /* 1 ! + \Z I iP\ < C'^e(l + t)- 1+2CN£ , |J|<JV-3 

that we will now prove. We will divide the solution of Einstein's equations into a linear and a nonlinear 
part with vanishing initial data h„ v = + w^ u . The estimate for the linear part follows from: 



(16.2) w 



V 



5* 



Lemma 16.1. If v is the solution of 

□ v = 0, v\ t=Q = vq, d t v\ t=Q = vi 

then for any 7 > 0; 

(16.5) (1 + t)\v(t,x)\ < Csup ((1 + \x\) 2+ ^(\ Vl {x)\ + \dv (x)\) + (1 + |x|) 1+ > (x)|) 

x 

Proof. The proof is an immediate consequence of the Kirchhoff 's formula 

v(t,x) = t / (yi(x + toj) + (v' (x + toj), w)) dS(ui) + / vq(x + too) dS(uj), 

J\tj\=l J \uj\ = 1 

where dS(u) is the normalized surface measure on S 2 . Suppose that x = re±, where ei = (1,0,0). 

Then for k = 1 , 2 we must estimate 

(16.6) 

dS(u) _ f 1 Cdu! ^ f 2 Cds 



1 + \ rei + ^|fc +7 1+ ^ r _ tui) 2 + f 2 {1 _ w 2)) C*+T)/2 - 7 1 + (( r _ £ + ta )2 + t 2 fi ) (fc+7)/2 

If A; = 2 we make the change of variables t 2 s = r to get an integral bounded by Ct~ 2 and if k = 1, we 
make the change of variables is = r to get an integral bounded by i -1 . □ 

To estimate the nonlinear part we use Hormander's L 1 — L°° estimates for the fundamental solution 
of □, see [HillLTj: 

Proposition 16.2. If v be the solution of 

Uw = g, w \ t=Q = dtw\ t=0 = 

then 

l^(«,y)l 

- /,. /.,,. 

\i\<. 



;i6.7) \w(t,x)\(l + t + \x\)<C^2 f [ \ Z * 9 ( S ' V )\ dyds, 

fc 2 Jo Jb?i + 8 + \v\ 



Let /i^ = ^ n« + «V where, 



(16.8) Uw piU = -H al3 d a dph ljLU + F tJiU (h)(dh,dh), «vlt=o = ^"V^o = °' 

and 

(!6-9) Q(v = 0, tvlt=o = /i^| t=0 , 9t^| t=0 = d t h lu\ t=0 - 



We have 



16.10) |Z 7 F^(/i)(a/i,^)| < C |<9Z J /i| +C7 ^ \^3-\QZ K h\, 

\J\ + \K\<\I\ \A + \K\<\I\ + ^ 
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s) 2Cn£ . 



and since H a P = -h af3 + 0(h 2 ), 

(16.11) |z J (if a ^ a ^v)l <c YVT]\ dzKh \- 

Now 

(16.12) ! \dZ J h\\dZ K h\{s,y)dy< ]T ||3Z J fc(s, Oll£ 2 < C(l + 

|/j<jV 

We write h = h° + h 1 and estimate 

f 16 13) / l^°( g ^)| 2 dv < A f* f° ^ < CM 2 

and by Corollary 113.31 

(16.14) t ^ {S ^l dy<C t \dh\s,y)\ 2 w( q )dy<C 2 N E 2 (l + t) 2C »* 



, U + kl) 2 

where w is as in (|16.2|) . Hence 

(16.15) J ^^\dZ K h\(s,y)dy<Ce 2 (l + t) 2C ^ 
It now follows from Proposition 116.21 that 

ft £ 2 j 

(16.16) \w^(t,x)\(l + t + \x\)< J {1 + s)1 _ 2Cn£ <Ce(l + s) 2CN£ 
which proves ()16.4p . 
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